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ABSTRACT 

< 

^ We report new photometry (R, JHK S ) and spectroscopy (optical, near-infrared) 

taken from 2009-2012 of the highly variable young stellar object PTF lOnvg (also 

p/j known as IRAS 20496+4354 and V2492 Cyg), and discuss the limited historical images 

of this source. We also report longer wavelength single-epoch photometry from WISE 

r£3 (mid-infrared) and CARMA (millimeter) that were taken during the 2010 outburst. 

^ Following the 2010 maximum of RpTF<13.5 m and subsequent fade, during 2011 and 

2012 the source underwent additional episodes of brightening to i?ppp~ 16 m , followed 
by several magnitude dimming events including extended faint states at RpTF>19.5 m . 

1 The infrared color-color and color-magnitude behavior strongly suggests that extinc- 
tion variations in the circumstellar environment - with AAy up to 30 mag measured, 

j> and exceeding 70 mag plausibly inferred from historical data - are responsible for the 

large-amplitude photometric variations. Further, a ~221 day period is derived from the 
infrared lightcurves; 2010 optical data taken when the source was near maximum bright- 
ness do not phase well to this period, however. Spectral evolution in the red optical and 
near-infrared includes changes in spectral slope and correlated changes in the promi- 
nence of TiO/VO (optical) and CO (infrared) emission that was seen in the early stages 
of the outburst. Changes in these features as well as in atomic permitted emission line 
strengths are anticorrelated with changes in forbidden emission line strengths which, 
^ along with H2 emission, dominate the spectra at faint epochs. Notably, night-to-night 

variations in several forbidden line doublet strengths and ratios are observed. We also 
present the first high-dispersion optical (R-, I-band region) and near-infrared (Y-band) 
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data for this source, in a variety of photometric states. We show emission line profiles 
of atomic species (H I, He I, O I, Ca II, Fe I, and Fe II) likely formed in an accretion flow 
and/or impact. The origin of zero- velocity atomic Li I A 6707 in emission is unknown. 
The forbidden lines ([N I],([N II] [O I], [S II], [Ca II], [C I], and [Fe II]) are all blueshifted 
and likely arise from an outflow /jet. Several of these lines are also seen in a shocked re- 
gion that is spatially offset from the continuum source position. Blueshifted absorption 
lines of the Na I D doublet, K I A 7665,7669 doublet, and the O I 7774 triplet, as well 
as blueshifted absorption components seen against the broad Ha and Ca II triplet emis- 
sion lines, similarly are formed in the outflow. The CARMA maps resolve a spatially 
extended outflow in mm-wavelength CO. We attribute the observed photometric and 
spectroscopic behavior to rotating circumstellar disk material located at approximately 
0.5-0.7 AU from the continuum source, causing the semi-periodic dimming. Occulta- 
tion of the central star as well as the bright inner disk and the accretion/outflow zones 
renders shocked gas in the inner part of the jet amenable to observation at the faint 
epochs. We discuss PTF lOnvg as a source exhibiting both accretion-driven (perhaps 
analogous to V1647 Ori) and extinction-driven (perhaps analogous to UX Ori or GM 
Cep) phenomena. 

Subject headings: stars: formation - stars: pre-main sequence - stars: individual (V2492 
Cyg; PTF lOnvg; IRAS 20496+4354) 



1. Introduction 



PTF lOnvg, in the North America Nebula region of recent star formation, was identified by 
the Palomar Transient Factory (PTF; |Law et al.||2009| |Rau et al.||2009[ ) on 2010 July 8, as an 



optical transient based on automatic discovery and classification codes (Bloom et al. 2011). The 



source was rapidly followed up by the PTF collaboration. We announced our findings concerning 



the 2010 outburst of this Class I type young star in Covey et al. (2011) where we presented optical 



and infrared light curves and multi-epoch optical and infrared spectroscopy. Therein, we made 
the analogy of PTF lOnvg (now called V2492 Cyg) to the behavior of V1647 Ori, an embedded 
Class I-type young star that was observed over the past decade to undergo several large amplitude 
photometric events on few year time scales. 

V1647 Ori and PTF lOnvg both displayed large-amplitude (~4-6 mag) photometric rises from 
their faint states that were similar in the early stages to the outbursts of FU Ori stars. However 
these particular objects and others like them do not have the spectroscopic characteristics of FU 
Ori stars; rather than being absorption line dominated (especially at high-dispersion), they are 
emission-line dominated with absorption seen in only a handful of blueshifted features that are 
diagnostic of strong winds. Further, the V1647 Ori-type objects do not remain in the elevated 
photometric state for the long time scale associated with FU Ori outbursts (estimated at roughly 
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a century); instead, their brightening episodes last only a few months to a few years and are 
characterized by large amplitude fluctuations on month to few month time scales. Members of this 
category possibly undergo repeated episodes of their outbursting (and/or extinction dominated) 
behavior at several year intervals. In this regard they are similar to the lower amplitude (~2-4 
mag) but repeating outbursts of EX Lup-type systems, which last months to >1 year each and 
repeat on few year to decade intervals. 

Whether small-scale, low-amplitude events such as EX Lup-type outbursts, larger scale, high- 
amplitude events such as the V1647 Ori-type events, or similarly large amplitude but also longer 
duration FU Ori-type events, the star/disk system in such outbursts is interpreted as undergoing 
an episode of enhanced mass accretion and associated mass outflow. The accretion mechanisms for 
the different categories of objects are likely related, and are attributed to instabilities in the inner 
accretion disks possibly associated with cyclic magnetosphere-disk interactions. For PTF lOnvg in 
particular, Covey et al. (2011 ) estimated a mass accretion rate of 2.5 x lO -7 M0yr -1 in the elevated 



state of 201CQ similar to the value estimated for EX Lup during its 2008 outburst (e.g. ( Juhasz et 



al. 2012)). 



Time- variable accretion, however, may not be the entire explanation for many large-amplitude 
young star variables. As demonstrated herein, time-variable extinction is clearly an important part 
of the PTF lOnvg interpretation, and may also play a significant role in the observed photometric 
behavior of many of the so-called "outburst" light curves highlighted to date in the literature. For 



example, GM Cep was discussed by Sicilia-Aguilar et al. (2008) as an EX Lup-type object, but 



later assessed by Xiao et al. (2010), Semkov & Peneva (2012), and Chen et al. (2012) as having 



extinction-dominated rather than accretion-dominated time series behavior. Both phenomena are 
perhaps simultaneously relevant, as we argue here for PTF lOnvg. Other well-known low-amplitude 
young star variables such as UX Ori, RR Tau, VV Ser, AA Tau as well as the "dippers" discussed by 



Morales-Calderon et al. (2011) and Cody & Hillenbrand (2010, 2011) also appear to be undergoing 
short time scale extinction events. Larger amplitude and - importantly - long period examples 
include KH15D ( [Hamilton etH|2012l >, WL4 QPlavchan et"alp058| >, and YLW16A QPlavchan et~ 



2010) which have found explanation in binary interactions with circumstellar disk material. 



The 2010 brightening of PTF lOnvg was independently detected by K. Itagaki and reported 
in 2010, August in the subscription service of the Central Bureau for Astronomical Telegrams as 



CBET 2426. Additional papers to date discussing the post-outburst source include those by Kospal 



et al. (2011) and Aspin (2011) who also present multi-color photometric, spectroscopic, and SED 



analysis of the source, often using the Covey et al. (2011) data. 



In this paper we present new data gathered by us since the publication of Covey et al. (2011). 



Intensive time series photometry shows that the source has continued its large amplitude and 
color photometric fluctuations. Time series spectroscopy demonstrates 1) continuum changes; 2) 



1 The accretion rates quoted in 



Aspin 



(2011 1 are a factor of 2.5-10 higher. 
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variation in the broad TiO / VO optical band emission that was detected for the first time during 
the 2010 outburst of PTF lOnvg, along with variation in other molecules such as CO and H2; and 3) 
drastic changes in the atomic emission lines indicative of accretion and outflow that are correlated 
(permitted lines) and anti-correlated (forbidden lines) with the photometric brightness. We also 
present high-dispersion spectral data; the line profiles are used to quantify the velocities relevant 
to the inflowing and outflowing material. An updated post-outburst spectral energy distribution 
is discussed, demonstrating that the source brightened relative to its historical spectral energy 
distribution in the mid-infrared, as well as in the near-infrared and optical, during 2010. A spatially 
extended, low velocity molecular outflow is detected. New high spatial resolution near-infrared 
imaging rules out the presence of stellar companions within several hundred AU. 

In ^2] we describe our new observations from 2009-2012. We then present our analysis of the 
multi- wavelength light curves (^3]), changes in the overall spectral energy distribution (Q, and 
analysis of the continuum, absorption line, and emission line spectroscopy including discussion of 
spatially offset forbidden line emission (^5]). In ^6] we present constraints on the source multiplicity 
from high angular resolution direct imaging. In ^7] we interpret our results and discuss the broader 
implications of our findings for accretion and extinction evolution of young stars. Finally in S|8] we 
summarize and conclude. 



2. Observations, Data Reduction, and Results 

Our photometric and spectroscopic monitoring of PTF lOnvg (V2492 Cyg) between 2009 and 
2012 includes continued monitoring in the R-band as well as J-, H-, and K s -band time series 
photometry. Almost thirty epochs of low dispersion spectroscopy at optical or infrared wavelengths 
were obtained. High dispersion spectroscopy in the red optical or in the 1 /im atmospheric window 
were acquired at eight epochs. In this section we describe the details of these various observations 
and resulting data sets. 



2.1. New Photometry 

2.1.1. Optical data from the Palomar Transient Factory 

Continued optical monitoring of PTF lOnvg was conducted as part of the regular observations 
of the greater North America Nebula field with the main PTF Survey Camera on the 48 inch Samuel 
Oschin Telescope at Palomar Observatory (hereafter P48). All observations were taken with 60 sec 



exposures in the -Rptf filter, a Mould R filter which is similar to the SDSS-r band (see Law et al. 



band (see 


Law et al. 


Covey et al. 


(2011 


); 



2009). P48 observations from 2009 August - 2010 December are presentet 
the new observations reported here started in 2011 February and continued through 2012 July. 
The typical cadence during the 2011 observing season was two observations separated by roughly 
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one hour obtained every night during astronomical bright time, though there was a period of high 
cadence observations during 2011 July with as many as eight observations taken in a single night. 
On 2012, April 18 a cadence of 3 observations was adopted. 

Brightness variations of PTF lOnvg are measured via a point-spread function (PSF) fitting 



routine as described in Sullivan et al. (2011). Briefly, in each image frame the PSF is determined 



from several field stars and the average PSF is then fit to the position of PTF lOnvg with each 
pixel weighted according to photon statistics. Saturated pixels are masked from the fit, enabling 
photometry on even moderately saturated point sources (relevant near the absolute peak in 2010). 
We note that the flux from PTF lOnvg may be slightly underestimated, particularly when the 
source is faint, due to oversubtraction of the complex background, which includes a bright arc of 
nebular emission ~7"to the east of PTF lOnvg. The flux calibration is performed relative to SDSS, 



based on nightly observations of SDSS calibration fields (see Ofek et al. 2012 for further details), 
and all reported magnitudes are given for the native -Rptf filter in the Vega system^] Finally, 
we note that while the photometric measurements presented here are consistent with those we 



previously published in Covey et al. (2011), the improved calibration and PSF model used here 



supersede our previous photometry. The optical light curve for PTF lOnvg is shown in Figure [T] 
and photometric measurements are reported in Table [TJ From 2011 June 08 through 2011 August 
21, and 2011 December 20 to 2012 April 28, PTF lOnvg was below the 3-a detection limit of the 
P48 images. In order to constrain the brightness of PTF lOnvg during these extended faint periods, 
we astrometrically aligned the P48 images and dynamically stacked them to determine the average 
brightness of PTF lOnvg over roughly two week intervals. 



2.1.2. Near-Infrared Data from PAIRITEL 



Continued near-infrared monitoring of PTF lOnvg was conducted with the 1.3 m Peters Auto- 
mated Infrared Imaging Telescope (PAIRITEL; |Bloom et al.|2006[ ) on Mt. Hopkins, AZ. PAIRITEL 
is a roboticized system using the former 2MASS southern hemisphere survey camera that employs 
two dichroics to observe simultaneously in the J, H, and K s bands. Observations were scheduled 
and executed via a robotic system. PAIRITEL is operated in a fixed observing mode in which 7.8 s 
double-correlated images are created from the difference of a 7.851 s and a 51 ms integration taken 
in rapid succession (see Blake et al.|[2008 ). The standard observing procedure involves taking three 
image pairs prior to dithering the telescope. 

The raw data from these images are reduced using standard IR reduction methods via PAIRITEL 
PIPELINE III and the flux for all sources is measured via aperture photometry using SExtractor 



(Bertin & Arnouts 1996). The final absolute calibration is determined using high signal-to- noise 



2 We convert from the AB mag system to the Vega system using the following transformation: m/{(Vega) = 
roji(AB) - 0.19. 



- 6 - 




50 75 100 300 400 500 600 700 800 900 1000 1100 



JD - 2455000 (day) 

Fig. 1. — Multi- wavelength lightcurve of PTF lOnvg with UT dates indicated above the figure. 
From bottom to top the data streams represent variability in the R-band (red; data from PTF) 
and in the J-, H-, and K-bands (black, purple, and blue respectively; data from PAIRITEL). 
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detections from 2MASS. When PTF lOnvg is very bright, corresponding to roughly m < 10 mag, 
it saturates the 7.851 s frames. This occurs in most of the K s band imaging and the brightest H 
band epoch. PIPELINE III produces "short-frame" mosaics consisting of reduced, stacked 51 ms 
images which can be used to measure the flux from bright sources (see also Bloom et al. 2009). 
The "short-frame" mosaics contain a few dozen bright 2MASS stars which we use to calibrate the 
photometric measurements of PTF lOnvg in the 2MASS photometric system. PAIRITEL has a 
known systematic uncertainty of a few percent on all flux measurements in each of the J, H, and 
K s bands (see Blake et al. 2008, Perley et al. 2010), which is larger than the typical statistical 
uncertainties. Following the method in Perley et al. we estimate the systematic uncertainties by 
measuring the scatter in the calibrated flux measurements of several (11 in this case) bright stars. 
Following this procedure we estimate systematic uncertainties of 0.025, 0.03, and 0.06 mag in the 
J, H, and K s bands, respectively. These systematic uncertainties are added in quadrature with 
the statistical uncertainties to provide the final uncertainties on each near-infrared measurement 
produced by PAIRITEL. We note that occasionally large thermal backgrounds prevent the detec- 
tion of sources in the K s band images produced by PAIRITEL, which has affected some of the 
epochs included in our analysis. The near-infrared light curve of PTF lOnvg is shown in Figure [T] 
and the photometric measurements are summarized in Table [2j which marks measurements made 
on short-frame mosaics as well as those affected by thermal background. 



2.1.3. Mid- Infrared data from WISE and Herschel 

The Wide-field Infrared Survey Explorer (WISE) observed PTF lOnvg as part of routine survey 
operations, obtaining several tens of scans in late May and early June of 2010. This was around the 
time of the first peak but before the absolute peak brightness in the optical light curve shown in 
Figure [l| The reported WISE photometry is partially saturated in bands Wl (3.35 ^m), W2 (4.60 
/im), and W3 (11.56 /im) with 12-15% of the pixels in the point-spread function in the saturated 
regime. The profile-fitted magnitudes account for the saturation while the aperture magnitudes 
are saturated and therefore lower limits to the fluxes. The W4 (22.24 /im) band photometry is 
unsaturated. 

Additional observations with WISE were obtained during the post-cryogenic phase of the 
mission, in late November of 2010. These would consist of only Wl and W2 data, but as the 
observations took place near the source's photometric minimum for the 2010 season they could 
be unsaturated. These data are not yet available in the WISE archive. Additional mid-infrared 
photometric epochs at similar wavelengths are being obtained with Spitzer in a sporadic fashion 
during 2011 and 2012 (program PI Abraham). 

Herschel has also observed PTF lOnvg in a program by A. Kospal and imaging photometry 
at 70, 100, and 160 /im (taken 2011 October 29, 2011 November 29, 2012 January 5, and 2012 
January 11) and 250, 350, and 500 /im (taken 2012 January 3) should be available soon. 
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2.1.4- Millimeter- Wave data from CARMA 

PTF lOnvg was observed in the 2.7 mm continuum, 12 CO J=l-0, and 13 CO J=l-0 with the 
Combined Array for Research in Millimeter- wave Astronomy (CARMA), located at an altitude of 
7200 feet in the Inyo mountains of eastern California. The observations used six 10.4 m diameter 
antennas and nine 6.1 m diameter antennas. Observations in the CARMA C-configuration, which 
provides baseline lengths between 26 and 370 m, were acquired on 2010 November 1 and 2 UT, 
near the 2010 photometric minimum. Observations in the CARMA D-configuration, with baseline 
lengths between 11 and 148 m, were obtained on 2012, June 22 and 23 UT. 

For the C-configuration observations, the local oscillator was set to a frequency of 107.7 GHz. 
The CARMA correlator contains 8 bands. Six bands were configured to have 487 MHz bandwidth 
to observe the continuum. The remaining two bands were set to observe the J=l-0 rotational 
transitions of 12 CO (115.271204 GHz) and 13 CO (110.201353GHz) with 0.98 kHz resolution (~ 
0.25 km s _1 ). Time-variable atmospheric and instrumental gains were calibrated by observing the 
source J2038+513 every 20 min. Between gain calibrator observations, PTF lOnvg was observed 
for 7 min, followed by a 7 minute observation of a second embedded star in the North American 
Nebula. The passband amplitude and phases were calibrated using observations of 3C454.3 on 
each night. The absolute flux was calibrated by observing Uranus on November 1 and Neptune on 
November 2; the estimated absolute flux uncertainty is 10%. The visibility data were calibrated 
using the MIRIAD reduction package. 

The D-configuration observations were similar to the C-configuration observations, with the fol- 
lowing differences. Two of the correlator bands were centered on CN, N=l-0 J=3/2-l/2 (113.490982 
GHz) and C 18 0, J=l-0 (109.78216 GHz), but these data are not presented here. The passband was 
calibrated with observations of 3C84. Flux calibration was performed by observing MWC 349 on 
the first night and Neptune on the second night. Finally, on the second night in D-configuration, 
a seven point mosaic was obtained to image the larger scale molecular cloud. 

Continuum images were generated by inverting the complex visibilities using "natural" weight- 
ing for optimal point-source sensitivity. The images were then deconvolved using the MIRIAD task 
mossdi. The RMS noise in the image measured in emission-free regions near the vicinity of PTF 
lOnvg is 0.18 mJy/beam, with a full-width-at-half-maximum (FWHM) beam size of 3.6" x 3.1". 
Figure [3] shows the 2.7 mm continuum map overlaid on a PTF optical image. A compact source 
with a signal-to-noise of 21 is detected in the field at a right ascension and declination of 20:51:26.2, 
+44:05:23.8 J2000, corresponding to the optical/near-infrared position of PTF lOnvg. A secondary 
compact peak in the continuum image is located 11" to the southwest; this peak has no counterpart 
in the optical image shown, or in archival ground-based or Spitzer images between 1.2 and 70 //m. 

The integrated 2.7 mm continuum flux toward PTF lOnvg is 5.60 ± 0.26 mJy (random error), 
obtained by fitting an elliptical gaussian to the image. Formally the source is extended relative to 
the synthesized beam with a deconvolved size of 2.7" x 1.5" and a position angle of 66°. However, 
when imaged with the C-configuration data only, the source is unresolved with a synthesized beam 
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Fig. 2. — Wide-field view showing velocity-integrated (3-6.25 km s _1 ) 12 CO (top panel) and 13 CO 
(bottom panel) images in the color scale compared with PTF R-band intensity in the contours. 
The location of PTF lOnvg is marked by the "x" towards the upper left. A remarkable alignment is 
found between the arc of Ha nebulosity (see also Figure [5]) and the eastern edge of the CO emission. 
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Fig. 3. — Left panel: CARMA 2.7 mm continuum image (contours) overlaid on a PTF R-band 
image (color scale) from 2010. Note the reversal in color vs contour representation compared to 
Figure [2j Solid contours begin at 3a with 3a increments, where a = 0.18 mJy/beam and the 
beam size is 3.6" x 3.1". Dashed contours start at —3a. The strongest 2.7 mm continuum source 
is coincident with PTF lOnvg. A second continuum source is located 11" to the southwest, but 
has no counterpart in optical or infrared images out to 70 fim. Right panel: Contour map of the 
12 CO J=l-0 integrated intensity between 6.5 and 11km s _1 overlaid on the PTF optical image. 
Contours begin at 3a with 3a increments, where a = 96 mJy beam -1 km s _1 and the beam size is 
1.7" x 1.4". The 12 CO emission at these velocities traces a red-shifted lobe of an outflow with PTF 
lOnvg at its apex. The corresponding blue shifted lobe of the outflow was not detected. 
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of 1.8" x 1.5" and an integrated flux of 3.2 ±0.30 mJy. Similar combinations of extended and point 
source components have been observed in protostellar sources (e.g., Jorgensen et al. 2004), and 
have been interpreted to represent a resolved envelope and a compact disk. The presence of these 
structures in PTF lOnvg are consistent with the Class I classification from the shape of the spectral 
energy distribution. 

The 12 CO and 13 CO spectral line images were generated using a Brigg's robust factor of zero to 
yield a beam size of 3.2" x 2.8". Since the emission is extended over several arcminutes, we used the 
maximum entropy program in MIRIAD (mosmem) to deconvolve the images. However, we caution 
that since no zero spacing data are available, the extended features are not reliably recovered in 
these images. Figure [2] shows images of the integrated 12 CO and 13 CO intensity between velocities 
of 3 and 6.25 km s _1 , which covers the main part of the cloud. 

PTF lOnvg, denoted by the cross in the figure, is coincident with a local maximum in the 
12 CO and 13 CO emission. Emission from the larger scale extends to the west and southwest over 
a spatial extent of at least 0.18 pc. The full extent of the local cloud has not been encompassed 
with the CARMA mosaic, however. The detailed spatial structure of the extended cloud differs in 
the 12 CO and 13 CO images, which we attribute both to differences in optical depth and the lack of 
sensitivity to extended emission with an interferometer. 

To the east, both the 12 CO and 13 CO emissions falls off sharply, suggesting that this structural 
feature has been robustly imaged with CARMA; PTF lOnvg is located along this ridge. The white 
contours in Figure 2 show the PTF optical image, which encompasses the Ha line. The Ha mission 
peaks to the east of the 12 CO and 13 CO emission, and closely follows the edge of the molecular cloud. 
Thus the eastern edge of the cloud appears to be illuminated by stars that are photo-dissociating 
the H2 on the cloud surface. 
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Fig. 4. — 12 C0 (solid) and 13 C0 (dashed) J=l-0 spectra in the vicinity of PTF lOnvg. Panels are 
labeled with positional offsets in arcseconds relative to the position of PTF lOnvg. Both the 12 CO 
and 13 CO emission contain strong emission at the cloud systemic velocity that peaks at the location 
of PTF lOnvg (see Figure [2]). The 12 CO emission also displays reshifted emission that extends to 
the south of PTF lOnvg (see right panel of Figure pi) that represents a redshifted outflow. 
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Figure Q shows the 12 CO and 13 CO spectra within ±2" of PTF lOnvg. Both 12 CO and 13 CO 
contain strong emission at a velocity of ~ 5 kms -1 , which corresponds to the systematic velocity 
of the cloud. As shown in Figure [2| the molecular emission at these velocities contains a local 
maximum at the position of PTF lOnvg. The 12 CO line profiles in Figure [4] also display weaker 
emission that is redshifted up to 10-15 km s -1 from the systemic cloud velocity. The right panel 
in Figure [3] shows an image of the redshifted image constructed from the C-configuration data. 
The 12 CO emission at these velocities traces an arc-like structure extending to the south, with 
PTF lOnvg at the apex. The morphology and the redshifted velocities indicate an outflow from 
the young (proto)star. No evidence is seen in the 12 CO spectra for a corresponding blueshifted 
component of an extended flow. 



2.2. Adaptive Optics Imaging 

To assess the multiplicity of PTF lOnvg at wide orbital separations, we acquired high spatial 
resolution near-infrared images on 2010 September 25 UT using NIRC2 (K. Matthews) and the Keck 
II adaptive optics system at Mauna Kea. The source was sufficiently bright during its outburst 
phase to serve as its own natural guide star at the time of our observations. We obtained 18 
unocculted, dithered frames in the K s band totaling 135 sec of on-source integration time. The 
NIRC2 narrow camera mode provided fine spatial sampling (10 mas/pix) of the stellar point- 
spread function (FWHM~ 65 mas for these observations). Raw frames were processed by cleaning 
hot pixels, subtracting background noise, and aligning and coadding the results. As discussed in 
section [6j no companions are found in these observations. 



2.3. New Spectroscopy 



Tables 1 and 2 describe the optical and infrared spectroscopy obtained by our collaboration of 
PTF lOnvg since its photometric outburst in 2010. Low resolution optical spectroscopy generally 
covering the red optical spectral region (>6200 A) were acquired at over 20 epochs with the facility 
spectrographs at one of: the Keck I 10m (LRIS; |Oke et al.|1995[ [McCarthy et al.|1998[|Steidel et al. 
2004] ), Palomar 5m (DBSP; |Oke Gunn(l982| |, or Lick 3m (Kast; pyHler fc Stone|l993l > telescopes. 



High dispersion optical data were also obtained at Keck I using HIRES (Vogt et al. 1994) on 5 



occasions. In the infrared, low resolution data covering approximately 1-2.5 /im were obtained at 



over 10 epochs using facility spectrographs at: the 3m InfraRed Telescope Facility (SPeX; (Rayner 



et al.|[2003 )), Palomar Observatory 5m (TripleSpec; ( Herter et al.|2008) >), and 3.5 m Apache Point 
Observatory (TripleSpec; ( Wilson et al.|[2~0 04)). High dispersion infrared data in the Y-band (1 /im 
atmospheric window) were also obtained, at Keck II using NIRSPEC (McLean et al. 1998) at 3 
epochs. 



The data acquisition strategy included sensitivity to obtaining observations at the parallactic 
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angle, and to proper acquisition of flux calibration standards. Spectral data processing followed 
standard procedures for bias level correction, flat fielding, spectral extraction including background 
subtraction, and flux calibration for the low dispersion data. Of note is the extended emission from 
a nebular arc region about 7" away from PTF lOnvg (see Figure [5] for morphology) that must be 
accounted for when perfoming background subtraction on spectra. The absolute calibration of the 
low dispersion data is typically good to ±20% and slightly better under photometric conditions. 
The resulting signal-to- noise ratios in the optical spectra range from 15-45, but lower in the 2011, 
late June to early August time frame when the source was optically very faint. For the moderate 
resolution infrared spectra, the signal-to-noise was typically ~250 in the K-band and ~40 in the 
J-band though varied by factors of a few depending on the weather and source brightness, and 
color at any particular epoch. 



3. Photometric Analysis 
3.1. Optical and Near-Infrared Variability 



The light curves published in Covey et al. (2011 ) included data obtained through 2010 Decern 



ber 14 (JD 2455545). Updated multi- wavelength light curves appear in Figure [T] and show that the 
source continued its large amplitude fluctuations throughout the 2011 and 2012 observing seasons. 
The data do not completely sample the photometric behavior, notably during the period when the 
field was too close to the Sun from approximately late December to early February of each year or 
during survey-designed gaps in the Rptf observing, as well as when the source was optically faint 
and below the PTF detection limits. Despite the incomplete sampling, the photometric evolution 
of PTF lOnvg is obviously dominated by large amplitude and variations over several month time 
scales, with additional chaotic behavior on shorter time scales. 



As described in Covey et al. (2011), PTF lOnvg was not a recognized optical or near-infrared 



source prior to its 2010 identification as an outbursting (possibly transient) source. It is noted in 
Aspin (2011), however, that the source is reported at faint optical magnitudes in the USNO-B1 
catalog in 1979, November [^] as well as in a 1997, October 30 Ha image (their Figure 2); it is 
also recorded in UKIDSS infrared data from 2006, June. These detections led us to investigate the 



3 In the digitized scans of all POSS images including the position of PTF lOnvg, the source is detected only once 
in a red plate taken in 1990 (see Figure [5]). The B and R mag measurements reported in Aspin (2011) correspond to 
the USNO-B1.0 source 1340-0394684, which has a cataloged position that is ~3.4" away from the optically measured 
position of PTF lOnvg. We note that the USNO-B1.0 catalog reports an epoch of 1966.2 for this source, which is the 
result of considering the Bl and Rl plates from the early 1950's (on which the reported photometry is in fact likely 
related to the nebulosity rather than the point source), and the R2 plate from 1990 (on which there is a true source 
detection) with the final position given as the geometric mean of these three epochs after taking a potential proper 
motion term into account (S. Levine, private communication). However, we were unable to find a POSS image of the 
field taken in 1979. 
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Fig. 5. — Optical images with north upward and east to the left covering 96" of the PTF lOnvg 
vicinity, spanning many decades. Rows contain blue, red, and far-red data from top to bottom, 
while the columns contain different imaging surveys, as indicated above the figure. The individual 
panels correspond to different epochs. PTF lOnvg is located at the center of the magenta circle 
and is not visible at most epochs, though was apparent on the 1990 R band plate. The intersection 
of the red cross hairs indicates the position of the USNO-B source identified in Aspin (2011) as the 
likely quiescent counterpart to PTF lOnvg but no source is detected at this position in any of the 
POSS I plate scans. 
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historical variability of PTF lOnvg via examination of images taken over several decades prior to the 
2009 optical detection by PTF. In Figure [5] we show the vicinity of PTF lOnvg from several sources 
including the digitized plate scans of POSS-I and POSS-II, previously unpublished optical images 
from the KPNO-0.9m, and SDSS. The apparitions of the source in 1990, 1997, and 2006 combined 
with the large number of nondetections in other imaging data sets spanning approximately 50 
years (see Figure [5]) hint that this young stellar object has displayed long term, large amplitude 
variability. Notably, although the 1997, October Ha detection of the point source exhibited by 



Aspin (2011) is convincing, the source was not detected in BVRI or Ha images taken only a few 
months earlier (1997, August). Considering the body of evidence, the average faint-state magnitude 
in the red optical is likely below about 22 mag, based on the deepest digital data from the late 
1990's. The source clearly has been intermittently detectable, however. 

In August of 2009 the object was captured by PTF at around -Rptf ~ 17.5 mag in what may 
have been the decline from an earlier local peak, with subsequent fading over approximately a 1 
month time scale to below ~ 19.5 mag, and thus non-detection by PTF. Then in 2010, after a four 
magnitude rise in less than 200 days, PTF lOnvg exhibited a first optical peak in July of 2010, 
then dropped in brightness by more than a magnitude peaking again approximately 40 days later 



in August of 2010 at an even brighter magnitude. The light curve published in Covey et al. (2011) 
exhibited a third peak and fall during early December of 2010. 

During the 2011 and 2012 observing seasons, several subsequent peaks that were typically 
narrower in time and fainter than the initial two 2010 peaks have been observed. Notably, the 
declines from each peak were sharp, approximately 0.1 mag day -1 , while the rise times were about 
half as steepQ The total optical variability amplitude observed during the 2011 and 2012 seasons 
was more than 4.5 magnitudes, with PTF lOnvg as bright as -Rptf ~ 15 mag in early 2011 and 
repeatedly reaching -Rptf ~ 16 mag during the rest of 2011 and 2012, but also repeatedly fainter 
than the P48 60-second detection limit (~ 19.5 mag in this high background field). Stacked images 
from the extended faint periods (e.g. between late June through late August in 2011 and again 
from late December 2011 through mid-March 2012) indicate that PTF lOnvg was on average during 



these times -Rptf ~ 20.4 mag (see § 2.1.1). Near- infrared detections continued through the periods 
of optical dimming and better sample the source fading (e.g. in 2011, mid-July), and rising (e.g. 
2011, early September) while it remained undetected optically. By 2011, November the source had 
brightened to -Rptf ~ 15.7 mag before once again fading below the P48 60-sec detection limit 
towards the end of the calendar year. In early 2012, PTF lOnvg remained fainter than -Rptf ~ 
19.5 mag but underwent a brightening eposode beginning in 2012, late March with a subsequent 
dip and rise again in late 2012, June and continued brightening through early August. 



4 The individually measured slopes for events lasting longer than a week are: 0.14, 0.14, and 0.14 mag day -1 for 
the drops following the peaks that occurred on ~JD 2455678, 2455884, and 2456086, respectively. For the rise times, 
the measured slopes are ~0.05, 0.05, and 0.07 mag day -1 characterizing the events starting on JD 2455652, 2455820, 
and 2456058, respectively. 
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Fig. 6. — Color-color and color-magnitude diagrams with the JHK S lightcurve data exhibited in 
Figure [T] shown as black points. Cyan pentagons show two measurements from Aspin (2011) while 
the gold star shows the 2003 measurement from UKIDSS. For reference, the dashed line indicates 
the locus of classical T Tauri stars as measured by Meyer et al. (1997) using the color transform of 



Carpenter (2001 ). When the overall features of the dataset are considered, the time variation of the 



black and cyan points is consistent with the effects of extinction following an interstellar extinction 
law from Indebetouw et al. (2005), plotted in dash-dot lines. The point-to-point variations, however, 
may require additional thermal effects beyond the dominant extinction effect. 
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The infrared light curves (Figure [TJ generally follow the fading and peaking behavior seen in 
the optical, with no detectable time lag between the optical and infrared valleys and peaks seen 
within our coarse sampling. However, the amplitudes are different between the RJHK S bands, and 
indicate large color variations. PAIRITEL obtains data simultaneously in the J, H, and K s bands 
allowing for precise measurements of color change in the near-infrared. A color-color diagram (left 
side of Figure [6]) indicates that during its large amplitude brightness changes, the color evolution 
of PTF lOnvg does not deviate strongly from the behavior expected from variations in line-of-sight 
extinction. The observed scatter around the expected relationship for time- variable obscuration by 
standard dust grains is roughly consistent with the typical observational uncertainties. A color- 
magnitude diagram (right side of Figure [6]) reinforces the finding that the source becomes redder 
when fainter in a manner consistent with extinction^ 

Assuming a standard interstellar extinction law and that all of the observed near-infrared color 
variations are the result of changes in extinction along the line of sight, we find that the observed 
excursions in near-infrared color-color and color-magnitude diagrams correspond to variations in 
extinction AAy > 30 mag (see also discussion in § [5j0 related to Figure[9]). Extending the analysis, if 
we adopt the position that all of the photometric variation oberved for this source can be attributed 
to line-of-sight changes in extinction (rather than to an outburst), we can compare the 2MASS non- 
detection at K s -band to the peak brightness observed with PAIRITEL to infer AAy ~ 70 mag. 

With evidence for large amplitude systematic changes in the photometry, we sought to quantify 
the time scales associated with the brightness changes. Repeated brightening and fading suggests 
periodicity. However, some general challenges for interpreting the lightcurve of PTF lOnvg include: 

(i) the inhomogeneous sampling that fails to resolve all the peaks and valleys in the light curve, 

(ii) the limited time baseline of the observations that may not properly resolve long term trends, 
and (iii) the likely possibility that the system is dynamic, meaning that periods, amplitudes, and 
relative phases could be changing cycle to cycle. 

To search for periodic signals in the variable brightness of PTF lOnvg we use a generalized 
Lomb-Scargle periodogram (Lomb 1976 Scargle 1982 Zechmeister & Kiirster [2009 ) and analyze 
each photometric band separately (see Richards et al.|2011 for more details on our implementation). 



3 We note that the same conclusion is not reached in 



Aspin 



(2011 1 based on an analysis of the UKIDSS detection 



from 2006 and seven additional observations made between 2010 July-2010 September. We find that indeed the 2006 
detection by UKIDSS (J — H — 2.85 mag; H — K s — 3.4 mag) lies far from the ensemble of 2010 outburst through 
2012 points measured by PAIRITEL, as shown in Figure [6] However, the deviation argued for in Aspin (20111 of 
the observations taken during outburst from standard reddening laws appears to be the result of an assumption that 



the unreddened origin for PTF lOnvg lies along the "classical T Tauri locus." As shown in both Covey et al. (20111 



and Aspin (20111, however, PTF lOnvg was not a classical T Tauri star prior to outburst, but rather a Class I-type 



source. As Class I sources have both envelope and disk emission contributing to the near-infrared colors, there is no 
reason to expect that the origin of the reddening vector should lie on the classical T Tauri star locus. This realization 
does not, however, explain the discrepant colors reported by UKIDSS. One possibility is that in its faint state the 
source becomes more dominated by scattered light which could explain the turn blueward in the J-H color. However, 
the H-K discrepancy seems to require redder colors at this epoch not bluer. 
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The extended period of non-detections in the -Rptf band would lead to a biased periodogram since 
the framework for the Lomb-Scargle analysis does not account for upper limits. As a result, we 
focus on the three near-infrared bands. The form of the long-term secular changes in the light 
curve is difficult to identify over the relatively limited timescales that we probe. Thus, we perform 
our Lomb Scargle analysis by simultaneously fitting for the periodic and long-term linear trend in 
the data. The observations do not clearly justify the use of a higher order polynomial, yet we note 
that simultaneously fitting for either a quadradic or cubic long-term trend does not significantly 
alter the results described below. 

In Figure [7] we show the periodogram for each of the three near-infrared bands. The peri- 
odogram exhibits significant power on long time scales. In each of the J, H, and K s bands we find 
a large peak in the periodogram corresponding to a period P ~ 221 days. The power at the 1 day 
period alias, an artifact produced by the fact that we observe nightly, is strong, but weaker than 
the power inferred for the astrophysically significant period at 221 days. The peak is somewhat 
broad in frequency, perhaps indicating only semi-periodic behavior from cycle to cycle. 

To quantify the significance of the observed peak in the periodogram we randomly resample the 
PAIRITEL flux measurements while fixing the cadence to match our observations. Following this 
we repeat the Lomb-Scargle analysis described above and measure the peak value in the resulting 
periodogram, while excluding the 1, 1/2, 1/3, and 1/4 d overtones of the 1 d alias. This procedure 
is repeated 10,000 times, which allows us to determine a "false alarm probability" that the observed 
peak in the periodogram is the result of noise by examining the rank ordered values of the peak 
significance from the 10,000 trials. In Figure [7] we show the 0.1% false alarm probability, meaning 
that a peak in the periodogram greater than this value has <0.1% chance of being the result of 
random fluctuations. In each filter the peak is well above this value, which clearly shows that 
the measured periodicity is significant. This can visually be confirmed based on the phase- folded 
trend-removed light curves shown in the inset of Figure [7j We interpret these results as tentative 
evidence that the emission from PTF lOnvg varies periodically on timescales of roughly 221 d; 
however, we caution that the signal may be driven by the particular peaks and minima that we do 
observe, specifically the two deep minima observed near JD 2455520 and 2455760. We emphasize 
the importance that future observations will play in determining if the observed signal in the 
periodogram is real. With higher cadence observations taken over longer time baselines we will be 
able to confirm if the signal in the current data set is fundamental to the system, or a peculiarity 
of the observing window through which we have measured the source. 

We note that the Rptf photometry does not phase well with the period derived from the 
near-infrared (JHK S ) photometry, especially in the 2010 season. This suggests that while extinction 
clearly plays a prominent role in the variability, other processes such as variable accretion (which 
would affect optical photometry moreso than infrared photometry) could also be an important 
physical process. The correlation of enhanced optical TiO/VO band emission (requiring heating 
and then cooling of molecular material) and permitted line emission (requiring accretion) with the 
bright state that we discussion in §[5] supports this explanation. Further, we again note concerning 
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Fig. 7. — Analysis of time scales inferred from the photometric time series. Lomb-Scargle peri- 
odogram for each of the JHK S passbands. A period of around 221 days, though with a broad 
frequency peak, is apparent. It is well above the signal level of the calculated 0.1% false alarm 
probability signal, shown as the horizontal dashed lines and calculated using 10 4 random re-shuffled 
instances of the data. Insets show the phased data, with the simultaneously fit linear trend removed 
for clarity of presentation. The colors indicate sequential observing seasons or other gaps in the 
data as follows: cyan = 2010, red = first part of 2011, blue = second part of 2011, black = 2012. 
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the historical imaging of this source, that while it was not visible to the limits of the 2MASS survey, 
not once over the past three years of our photometric monitoring has it been fainter than K s ~ 
11 mag. As discussed above in the context of possible variations in extinction of Ay ~ 70 mag, 
the recently bright K s -band magnitude could alternately support the importance of an outburst 
component, despite variable extinction clearly being part of the picture. 

We can also consider the derived periodicity in the context of the intermittent visibility of 
PTF lOnvg in historical images (Figure [5]). Projecting backwards in time the 221 day period 
derived here from data taken over the past few years, we find that there should have been peaks in 
the vicinity of both the 1990, September DSS and 1997, October UKIDSS epochs when the source 
was apparent (see discussion above) . Similarly, the source would have been faint during the KPNO 
0.9m and the SDSS image epochs shown in Figure [5] We caution, however, that even a few day 
error in the period would lead to dramatically different conclusions when projected back over these 
several decades. 

Should the periodic signal be confirmed, that in conjunction with the evidence that the pho- 
tometric variations are largely due to changes in extinction along the line of sight (Figure [6]) would 
naturally suggest that the recently observed long period variability arises in a dusty structure char- 
acterized by a dynamical time T^ yn ~221 day. If this is a rotating disk around a 0.5 or 1 Mq 
star, it would indicate higher opacity material in a clump or warp located around 0.5 or 0.7 AU 
from the star. The implied orbital velocity for a circular orbit (25 km/s) can be compared to the 
characteristic time scale for the photometric rises and declines that would correspond to repeated 
clearing and obscuration events. Refering to Figure [TJ the relevant rise and fall times are a few tens 
of days. This time scale relative to the observed period is similar to the ratio between the derived 
location and the orbital circumference (2ir), resulting in a size scale for the obscuring material that 
is thus similar to its distance from the star, namely ~0.5 AU. 



3.2. Mid-Infrared Variability 

The data acquisition strategy for the WISE mission included multiple scans of each point on 
the sky, which provides time series data on time scales that vary with source ecliptic latitude. For 
the scans covering PTF lOnvg at moderate declination, the time series data are rather concentrated, 
with several observations obtained within just a few several-day-long time frames. The observations 
took place around the time of the first peak and before the absolute peak brightness exhibited in the 
optical light curve (Figure [I]). PTF lOnvg is flagged in the WISE catalogs as a potential variable 
in the W3 and W4 bands based on 26 individual measurements. Examination of the photometry 
from these individual scans reveals nothing systematic about the behavior in time between e.g. 
the 2010, May 27 and 28 WISE data vs the 2010, June 4 WISE data. Scatter among the many 
observations taken so close in time that is larger than the quoted errors likely leads to the variable 
source flag; it is not clear whether the source is actually varying on sub-day time scales, however. 
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Comparing the WISE photometry to previously reported measurements at similar wavelengths 
from Spitzer, MSX, and IRAS, illustrates that the source brightened in 2010 at mid-infrared wave- 
lengths as well as at the near-infrared and optical wavelengths highlighted in Covey et al. (2011) 
and discussed in detail above. Specifically, the Spitzer fluxes from 2006/2007 reported by Rebull 
et al. (2010) at both 3.6 //m and 4.5 /im were a factor of more than ten lower than those recorded 
by WISE in 2010 at 3.35 /mi and 4.6 //m. The WISE 22 /im measurment can be compared with 
data from not only Spitzer but also from earlier MSX (21.3 //m) and IRAS (25 /xm) measurements. 
While there are wavelength and beam size differences between these various missions to consider, it 
is clear that the 2010 WISE photometry is brighter by a factor of ~4 compared to the 2006 Spitzer 
24 /im data, with the earlier missions all reporting photometry at intermediate flux levels. 



4. Spectral Energy Distribution Variability 



In Figure |8| the outburst spectral energy distribution presented by Covey et al. (2011) is 
updated to include the photometry from WISE and CARMA. 

Using the CARMA continuum flux measurement of 5.6 mJy at 2.7 mm, a total mass (dust 
+ gas assuming a ratio of 1:100) of 0.06 M is calculated. This assumes a distance of 520 pc, a 
dust opacity of 0.01 cm 2 g _1 at 107 GHz (Beckwith et al. 1990), that the continuum emission is 
optically thin and isothermal with a temperature of 20 K, and that free-free emission is a negligible 
contribution. This value is on the high end of circumstellar masses measured for young stars and 
consistent with the Class I-type SED of the source. 

The post-outburst CARMA measurements at 2.7 mm can be compared with the pre-outburst 
data from the Bolocam Galactic Plane Survey at 1.1 mm (Rosolowsky et al. 2009) to search 
for variability at millimeter wavelengths. The Bolocam survey reported a 2a detection at the 
position of PTF lOnvg with an integrated flux of 148 ± 74 mJy, although the integrated flux may 
be underestimated by 50% (Rosolowsky et al. 2009). The 1.1 mm detection is consistent with a 
point-source at the 33" angular resolution of Bolocam. The higher resolution CARMA 2.7 mm 
continuum image resolves the continuum emission into two compact sources and a diffuse extended 
component (see Figure [3]) . The integrated flux of these three components measured with aperture 
photometry is 11.3 ± 1.4 mJy, where a 10% calibration uncertainty is assumed. While additional 
extended emission may be "resolved out" with CARMA, we nonetheless derive a spectral millimeter 
spectral slope (F u oc v a ) of a = 2.9 ±0.6. Assuming the millimeter continuum emission is optically 
thin, the expected slope is a ~ 4 for interstellar medium grains, and a ~ 3 if grains have grown 
to approximately millimeter in size. Thus within the considerable measurement uncertainty, the 
observed spectral slope is consistent with that observed in circumstellar disks and envelopes, and 
we find no evidence for variation in the millimeter flux toward PTF lOnvg. 



As demonstrated above, however, there has been significant variability at mid-infrared wave- 
lengths. The 2 — 25/mi spectral slope, originally defined for ground-based K-band and IRAS data 
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Fig. 8. — Spectral energy distributions updated from those presented in Covey et al. (2011) to 
include WISE and CARMA photometry. Bottom panel shows literature data from 1980-2009 
while upper panel shows outburst photometry from 2010. Optical and near-infrared photometry 
is selected to be closest in time to plotted mid-infrared and millimeter photometry. At WISE 
wavelengths, both the saturated and the profile-fit corrected magnitudes are shown. Downward 
facing triangles are upper limits and upward facing triangles are lower limits while squares and 
circles (large beam IRAS measurements in the left panel) are detections. As discussed in the 
text, the 2010 brightening of the source had mid-infrared as well as near-infrared and optical 
manifestations. 
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by Lada (1987) as a = d log \F\ / d logA, can be computed in the pre-outburst and post-outburst 
stages. This could inform us as to whether there have been differential brightness changes over this 
wavelength range indicative of changes in either the geometric structure or the radiative transfer 
processes in the circumstellar environment. Alternately, there could be relative brightness changes 
consistent with the variable extinction interpretation of the near-infrared photometry. 

Using the data shown in Figure [8] and with the caveat that there are not ideal matches in the 
timing between available K-band and mid-infrared photometry, we find indeed a dramatic change 
in the spectral slope of the source from pre-outburst to post-burst. Regardless of whether IRAS 
25 /j,m, MSX 21 /im, or Spitzer 24 ^irn photometry (taken over a time span of three decades) is 
compared to the 2MASS 2 upper limits from 2000, we find a spectral slope above 2 (in the 
range >2.1 to >2.5) whereas using the UKIDSS measurement in 2007 produces a spectral slope 
in the range 1.6-2.0, in either case indicating a "Class I" -type spectral energy distribution. The 
post-outburst measurements from PAIRITEL and WISE 22 /im produce a spectral slope ~0.3, still 
"Class I" but with a much flatter spectrum. This dramatic change in spectral slope is due to the 
fact that although there was demonstrable source brightening in the mid-infrared, by factors of 
at least several to ten, the near-infrared and optical brightening significantly exceeded that in the 
mid-infrared. 

Having established that the pre-outburst and post-outburst spectral slopes were different, 
though both indicative of a Class I-type SED, a key aspect of the outburst interpretation is to 
determine if the dramatic slope change is reflective only of the large changes in line-of-sight ex- 
tinction that were subsequently observed for this source, or if there was a true luminosity change. 
Unfortuntely we have no means of calculating the source luminosity independent of the extinction 
effects so can not address this issue. However, we note that a factor of four variation at 25 fim 
would be equivalent to a change of 1.5 mag in the extinction at this wavelength and correspond 
roughly to Ay of ~30 magnitudes. This is close to the maximum extinction change inferred from 
variations in the near-infrared colors (see Figures [6] and [9J . 



5. Spectroscopic Analysis 

Our new optical and near-infrared spectra at low dispersion can be considered along with the 



spectroscopic data published by Covey et al. (2011) that sampled the source both before and after 



its absolute brightness maximum in 2010, August. The 2011 and 2012 spectra sample epochs in the 
lightcurve of PTF lOnvg when the source was near local maxima in brightness as well as when it 
was deep in extended faint states e.g. during summer/fall of 2011 and fall/winter of 2012 (Figure [T]). 
The high-dispersion spectra also sample the lightcurve at relatively bright as well as much fainter 
photometric states. 

Figure [9] illustrates the ensemble of spectral epochs relative to the optical extinction variation 
over the time series, estimated by calculating the color variation from the bluest color observed (on 
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Fig. 9. — Differential optical extinction, A Ay, measured from the color changes in Figure [6] relative 
to the bluest photometry captured in our near-infrared time series and an assumed extinction law. 
Error bars are calculated assuming the typical photometric errors in Table [2j The alphabetical 
labels indicate the timing of our spectroscopic follow-up of the 2010 brighting of PTF lOnvg, with 
"H" representing high-dispersion optical data, "L" low dispersion optical data, "I" indicating low 
dispersion near-infrared data and "IH" high-dispersion near-infrared data, all as described in the 
text. We note in particular the spectroscopic coverage during the faint minimum when the central 
source was highly extincted. 
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JD 24555450) and applying the Indebetouw et al. (2005 ) extinction law. We note the high density of 



spectral coverage obtained during the faintest, reddest state of the PTF lOnvg light curve in 2011, 
summer. Considering the strongly detected photometric periodicity, we note that given the sparse 
spectral sampling we do not address directly the evidence for spectroscopic periodicity. However, 
we do observe a seemingly repeating pattern in the spectroscopic signatures that are discussed in 
detail below. Specifically, the forbidden lines are most prominent during the high extinction epochs 
represented in Figure [9| Further, spectra taken when PTF lOnvg was in one photometric cycle are 
very similar to those taken in earlier photometric cycles at similar phase; for example, repeatedly 
when brightening, the forbidden lines are clearly present but quite weak relative to their dominance 
in the faint photometric states. 

Our discussion of the spectral time series data is divided into presentation of the continuum 



behavior (illustrated for the optical in Figure 10 and for the infrared in Figures 11 and 12) in 



section 5.1, and then description of the lines. The line discussion is subdivided into wind-indicating 
absorption lines in section 5.2 and then accretion- and outflow- diagnosing emission lines in section 
5.3. The various emission species are further separated into molecular (section 5.3.1), permitted 
atomic (section 5.3.2), and forbidden atomic (section 5.3.3) lines. Evidence for spatially extended 
forbidden line emission is presented in section 5.3.4. A wide range of thermal and mechanical 
processes are probed with the array of available spectral diagnostics. The molecular lines measure 
gas having temperatures up to several thousand Kelvin. The atomic lines diagnose temperatures 
usually in the range ~5000 - 20,000 K (but up to 40,000 K), and densities n e « 10 4 - 10 9 cm" 3 
(forbidden lines) and nn ~ 10 9 - 10 12 cm -3 (permitted lines). Our main focus is on describing the 
evolution in time of the line presence, line ratios, and line profiles in the context of mass accretion 
and outflow. 



5.1. Evolution of the Optical and Infrared Spectral Continuum 



The large amplitude photometric changes exhibited by PTF lOnvg are mirrored by large 



changes in the spectral continuum. This is illustrated in Figure 10 showing the time series of 



optical spectroscopy within the wavelength range R and I-band regions and in Figure 11 showing 
the infrared spectra in the Y, J, H, and K-band atmospheric windows. 



Figure 10 shows the flux calibrated spectra illustrating the spectral continuum changes that 
accompany brightness changes of the source, and the normalized spectra illustrating the line-to- 
continuum changes. In this latter representation, lower signal-to-noise ratio spectra were taken 
during epochs when the source was fainter and show the presence of numerous prominent forbidden 
lines that are not seen in the brighter photometric states having stronger continuum emission. As 



discussed in Covey et al. (2011) the optical continuum shape of PTF lOnvg has varied, along with 



the photometry. In particular the TiO/VO emission features, which are broad in wavelength and 
thus act like a continuum, have varied in strength relative to the nearby out-of-band continuum 
regions. They also seem to have varied in shape, i.e. roundness vs flatness. We discuss these 



-27- 




6500 7000 7500 8000 8500 9000 

Wavelength [Angstroms] 




Wavelength [Angstroms] 

Fig. 10. — Top panel: Individual flux calibrated optical spectra. Brightness mimics the photometric 
light curve. Data from 2011, July 4 and July 5, and from 2012, January 1 are excluded due to 
significant noise. The spectrum from 2010, August 11 is not corrected for telluric absorption (the 
prominent atmospheric B and A bands have been smoothed over). Bottom panel: Normalized 
spectral time series plotted on a log scale (in order to emphasize the emission lines). Unmarked 
weaker emission lines are generally low excitation transitions of Fe I or Fe II. Note the slightly wider 
spectral range than in the top panel, highlightiing additional forbidden lines. Spectra are ordered 
by date with some dates having no accompanying spectra due to significant continuum noise. Both 
representations show that during fainter states (summer of 2011 and winter of 2012), the forbidden 
emission lines are more prominent. 
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changes in more detail in § [5}3.1 below when describing variations in the broader set of molecular 
emission features exhibited by PTF lOnvg. 

In the infrared, the continuum also exhibited differential brightening and fading. Along with 
the time series spectral observations shown in Figure 1 1 , we highlight the near-infrared spectra from 
the brightest and nearly faintest epochs in Figure 12 The variations in PTF lOnvg's near-infrared 
spectral continuum over the 2010-2011 seasons appear consistent with the potential extinction 
variations discussed in section 3.1 based on the PAIRITEL photometric monitoring data, and 
illustrated in Figures [6] and [9j To verify this assertion, we used an Ry=3.1 Fitzpatrick (1999) 
extinction law [j to de-redden each moderate resolution near-infrared spectrum according to the 
photometrically-derived Ay estimates in Figure [9| as interpolated onto each spectroscopic epoch. 
While the flux densities of the observed spectra differ by more than an order of magnitude, the 
flux densities of the de-reddened spectra agree to within a factor of ~2. Similarly, the spectral 
slopes of the de-reddened spectra also agree well generally, with the exceptions of the spectra 
acquired on 2010 November 26, 2011 June 26and 2011 July 15, when the source was deep in (two 
different) photometric minima. Considering the uncertainties imposed by the non-simultaneity of 
the spectra and the photometric data from which the extinction estimates are derived, this level 
of agreement suggests that extinction variations account for much of the near-infrared photometric 
and spectroscopic continuum variability within our time-series data. 



5.2. Evolution of the Wind-Indicative Absorption Lines 

Our spectra contain a number of lines sensitive to outflowing material. The early optical 



spectrum of PTF lOnvg shown by Covey et al. (2011) exhibited only the Na I D, K I A 7665 



7699, and O I A 7773 lines in absorption, all with blueshifted velocities relative to the expected 
line center wavelengths. All of these lines are either resonance lines or have a metastable lower 
level, so they readily go into absorption. Absorption was seen against the broad TiO/VO molecular 
emission continuum mentioned above, and there were no emission components to the profiles. The 
blueshifted absorption lines remained relatively constant with respect to a normalized continuum, 
i.e. in equivalent width, as the source varied photometrically over the 2010, 2011, and 2012 seasons. 
The equivalent widths of these lines hovered around 2-3 A over most of this time period, though 
possibly somewhat higher (3-4 A) towards the end of 2011. Between 2011 June 3 and 29, however, 
these absorptions disappeared as the continuum faded, perhaps only as they became undetectable 
in the lower signal-to-noise ratios characteristic of these spectra, or perhaps intrinsically if the wind 
strength decreased. By 2011, late August, when the source was again near a peak in brightness, the 



The 



Fitzpatrick 



( 1999 



reddening presecription (employed here to the spectra because it is smoothly varying 



with wavelength) and the Indebetouw et al. ( 2005 1 extinction (used earlier with the photometry since it is determined 



empirically using the same filters) are the same to within 5% in the near-infrared; this difference is roughly the same 
as the estimated error in the flux calibration of our near-infrared spectra. 
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Fig. 11. — Flux calibrated moderate resolution (R~2000-3000) near-infrared spectra of PTF lOnvg 
during 2010 and 2011. The majority of the spectra possess similar underlying continuua and 
emission line components. While the absolute flux densities of the observed spectra vary by a 
factor of ~10, once the spectra have been corrected with a standard Ry = 3.1 extinction law for 
the photometrically derived reddening estimates shown in Figure |9j the spread in flux reduces to 
a factor of 2, within expectations given the non-simultaneity of the spectra and the photometric 
extinction estimates. The spectra obtained on 2011 June 26 and July 15, when PTF lOnvg was in 
a particularly faint state ( J ~ 17.5 mag), lack the numerous atomic (e.g., H, Na, Ca, Mg, etc.) and 
molecular (e.g., CO, H2O) emission features that dominate the brighter spectra. The few emission 
features that remain visible in the faint state (e.g., H2, [Fe II]) likely arise in the jet/outflow, 
suggesting that the extinction may be sufficiently high to completely obscure the emission from 
near the photosphere, accretion region, and inner disk. 
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Fig. 12. — Comparison of PTF lOnvg near-infrared spectra during the brightest (2010 Septem- 
ber 23) and second faintest (2011 June 26, adopted because it is morphologically similar to the 
faintest 2011 July 15 spectrum but extends to bluer wavelengths) epochs captured in our moderate 
resolution near-infrared spectral monitoring. Numerous Fe I and other permitted atomic lines are 
prominent at ~1.15-1.2 microns in the brighter spectrum but absent in the fainter spectrum, which 
exhibits primarily H2 and [Fe II] . 
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Fig. 13. — Time series optical spectra in the same color scheme and order as Figure 10 showing Ha 
as well as the [N II], [O I], [S II], and [Fe II] forbidden line doublets. Lower signal-to-noise spectra 
have been smoothed. Given the large variations in line-to-continuum ratio over the time series, the 
ordinate is on a log scale. Note the variation in the doublet line ratios and the slightly blueshifted 
line centers which are apparent even at low resolution. 
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Fig. 14. — Time series spectra in the 7275-7474 A region showing the [Ca II] doublet as well as 
[Ni II], and [Fe II] lines that are prominent during faint photometric states. Lower signal-to-noise 
spectra have been smoothed, and neither of the 2011 July 4 or 5 spectra are shown. 
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line absorption from winds was again apparent at the same equivalent width level and is seen in all 
subsequent "bright state" spectra through mid-2012. The lack of change in absorption equivalent 
width is consistent with the extinction variation interpretation of the continuum behavior, implying 
that the large extinction variations arise exterior to the wind absorption region. 



Notably, the high-dispersion data from December of 2010 (Figure 15) shows that the Na I 
D doublet absorption was fully saturated down to zero flux, with a boxcar like profile extending 
blueward to ~ -275 km/s. The K I A 7665, 7699 doublet, similar to the He I A 10830 profile 
discussed below, has a very narrow absorption close to zero velocity yet slightly blueshiftecQ plus 
a separate multiple-component blueshifted absorption profile extending from approximately —50 
to —250 km/s. The depth is approximately 40% of the continuum level, shallower than the Na I 
D doublet absorption. The O I A 7773 triplet absorption profile is hard to decipher due to the 
component blending, but it appears approximately neutral in velocity and extends to appoximately 
60% of continuum depth; approximately 130 km/s of the line width is due to the triplet blending. 



We note that it is claimed in Aspin (2011) that the O I A 8446 line also exhibited a blueshifted 



absorption profile in 2010, September, but we see only emission from our low resolution data taken 



in August and September (Figure 10 ) and in our high resolution data taken in December (Figure 
15 ) of that year. 



The time evolution of the optical wind-indicative absorption lines can not be well-documented 
given that source faintness precluded the detection of continuum and thus any absorption against 
it in many of our spectra. In the HIRES data from 2011, December, however (taken after the 
source had undergone several episodes of photometric fading and brightening), the K I A 7665, 
7699 doublet and the Na I D doublet absorption profiles appear to have approximately the same 
form and terminal velocities as in the 2010, December spectrum. 



Edwards et al. (2006 ) have discussed the utility of the He I A 10830 line (also with a metastable 



lower level) as a wind indicator, demonstrating that its sub-continuum absorption over a range of 



velocities probes both disk winds and stellar/polar winds in young stars. Kwan Sz Fischer (2011) 
derived the physical conditions traced by He I (density njj ~ 10 9 cm -3 , temperature T ~ 10,000 K, 



excitation through UV photoionization) as well as other accretion/outflow tracing lines. Kurosawa 



& Romanova (2012) demonstrated in "observation" of MHD accretion/outflow radiative transfer 



models the sensitivity of this line to the high energy radiation field, namely Lx and Tx, and were 
able to produce line profiles having redshifted and blueshifted absorption and emission features 
consistent with observations of young accreting stars. The early infrared spectrum of PTF lOnvg 



shown in Covey et al. (2011) exhibited only He I A 10830 in absorption. The line was blueshifted 



relative to line center. In the left panel of Figure 17 we show the time series behavior of this line at 



low dispersion. While the range of resolutions of the spectra in Figure [17] do influence the details 
of the observed profile somewhat, the He I line appears to strengthen systematically over time, 



7 While this might be attributed to an interstellar contribution, the similar velocity of the narrow K I and narrow 
Na I to the non-interstellar narrow He I component is noteworthy. 
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Fig. 15. — Optical line profiles from 2010, December with vertical dotted lines indicating zero 
velocity. In each panel the black profile represents the line indicated first in the panel legend; 
blue/green profiles, if present, indicate the second/third line in the panel legend (for example, 
in the Ca II panel the black, blue, and green profiles are those of the AA8498, 8542, and 8662 
transitions). The Ca II triplet lines are broader than other atomic emission lines and likely have 
both accretion-related and outflow-related components; note that the red sides of these lines are 
each contaminated by H I Paschen lines, with the panel immediately below the Ca II panel showing 
another Paschen line on the same scale. Highly blueshifted emission is apparent in the forbidden 
emission lines of Fe II (see also Figure 16). Permitted Fe I, Fe II, and O I emission, by contrast, 
is centered at nearly zero velocity and relatively narrow, with the Fe I lines having the narrowest 
profiles. Blueshifted absorption is seen in the wind-sensitive absorption lines of K I and Na I with 
the blue sides of the Ca II triplet lines also exhibiting some (sub-continuum) absorption. The O I 
A7774 triplet is seen in absorption as well. Unusually, Li I A6707 appears to be in emission (see 
text). 
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Fig. 16. — Comparison of normalized optical Balmer and forbidden line profiles as identified above 
the top row. Spectra from bottom to top show data taken on 2010, December 13, 2011 June 15, 
2011 June 28, 2011 December 9, and 2012 January 6. The velocity range along the abscissa is 
— /+300 km/s. The ordinate scaling changes for the different dates (rows) but the relative scaling 
along each row is the same so that comparison of line ratios from row to row is meaningful; for 
example, O I to [Fe II] is high when the source is bright in 2010 December and 2011 December, 
but low in 2011 June and 2012 January when the source is faint. Black and blue lines indicate, 
respectively, the first and second wavelengths as labelled at top. Balmer line emission is visible over 
the same range of blueshifted velocities as the forbidden lines during faint photometric states, but 
there is an additional redshifted emission component during bright photometric states. Variation 
is seen in the ratio of the [O I] 6300 and 6363 lines and the [S II] 6717 and 6731 lines; for reference, 
the ratios are as expected for optically thin emission in the bottom row. 
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independent of the photometric state of PTF lOnvg. The line increases from ~3 A in equivalent 
width and a maximum depth of 60% of the continuum in 2010 July to ~8 A equivalent width and 
a maximum depth of ~30% of the continuum in 2011 October. 

These dramatic changes are also evident in the high-dispersion He I A 10830 profiles shown 
in Figure 18 The NIRSPEC spectra were obtained in 2010, July, when PTF lOnvg was near 
its absolute brightness maximum (which occurred in 2010, August) in 2011, September, when the 
source was in its prolonged optical minimum state, and in 2012, May as the source was rising towards 
a photometric peak. The change exhibited is consistent with the trend observed in the lower- 
resolution spectra. The faint-state 2011 spectrum is close to a classical P Cygni profile, with some 
sub-structure in each of the redshifted emission and the blueshifted absorption components, along 
with a narrow, slightly blueshifted absorption. This low velocity narrow absorption component to 
the He I A 10830 profile is present in all three spectra and may be the same as that seen in the 
K I profiles. In both brighter state spectra, the redshifted emission component at about +50 km/s 
that was present in the faint state is gone. However, the two brighter state spectra are different 
in their broad blueshifted absorption components. Relative to the 2010 bright state spectrum, the 
2011 faint state and 2012 bright state blueshifted absorptions both are deeper (25% of continuum 
vs 60% of continuum) and broader (about —100 to —300 km/s versus —50 to —250 km/s). 

In summary, the blueshifted absorption exhibited in the Na I D, K I A 7665, 7699, O I A 
7773 lines and He I A 10830 line, as well as blueshifted absorption components to the broad Ha 
and Ca II AA 8498, 8542, and 8662 emission profiles (discussed in § [5j3 below) are indicative of a 
strong outflow/wind. The absorption must occur in a region outside of the optical and near- infrared 
continuum in the case of the pure absorption features, and outside the region producing the broad 
Ha and Ca II line emission (mainly the accretion column but also the wind itself) . 



5.3. Evolution of the Accretion and Outflow Diagnosing Emission Lines 



At low resolution, Figures 13 and 14 show detail relative to Figure 10 of some of the optical 
emission line behavior of PTF lOnvg during 2010 and 2011. Among the rich optical emission line 
spectrum are permitted atomic species of: H I, O I, Fe I, Fe II, Ca I, and Ca II. Notably there is 
no evidence of optical He I. There are also (strong) forbidden emission lines in the optical from: 
[C I], [O I], [Fe II], [Ni II], [Cr II], [Ca II], and [S II]. Several of these species are rarely seen in young 
star spectra, even those which are emission-line dominated. Figure [17| illustrates the infrared time 
series of emission line behavior. In the infrared, we see permitted atomic species of: Al I, Ca I, 
Fe I, K I, Mg I, Na I, O I, Si I, and Ti I as well as forbidden [Fe II], [N I], and [S II]. 

In both the optical and infrared, there was notable weakening of permitted atomic line emis- 
sion when the continuum faded. The forbidden atomic line species, on the other hand, became 
significantly more prominent in the faint state spectra; this prominence is due not only to better 
contrast relative to the continuum as it faded, but also to variation (by up to a factor of two in 
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Fig. 17. — Time series of continuum normalized moderate resolution near-infrared spectra, isolating 
representative permitted and forbidden transitions. The two spectra with the highest photometric 
extinction estimates clearly lack emission lines that are thought to trace hot, dense circumstellar 
material (i.e., H I, CO, Mg I); by contrast, these spectra are dominated by lines thought to trace 
shocks (H2) or photo-excitation ([Fe II]) in lower-density circumstellar material. Considering lines 
with similar wavelengths but that trace different physical regimes (e.g., A2.12 H2 vs. A2.165 H I) 
systematic variations are measured as a function of the photometrically estimated extinction; the 
hot, dense tracer dominates when PTF lOnvg is relatively bright, and the lower-density tracer 
grows in prominence as PTF lOnvg fades. Measurements of the absolute line fluxes reveal that the 
lower-density tracers maintained a near-constant brightness while the lines tracing hotter, denser 
material faded along with PTF lOnvg's continuum. We interpret these aspects as evidence that 
the changes in line-of-sight extinction are taking place exterior to the protostellar photosphere and 
inner disk (and thus obscuring the continuum and hot, dense emission line regions) but inward of 
the lower-density regions emitting in H2 and [Fe II] . 
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Fig. 18. — Line profiles in the 1 /nm region as they appeared in 2010 July (black), 2011 September 
(red), and 2012 May (blue). Top row shows a He I and two H I lines. Bottom row shows metallic 
forbidden lines (two [N I] doublets near 1.0401 /jm and a [S II] line at 1.029 /im) and a permitted 
line (Si I 1.061 /mi). As in the high dispersion optical spectra, the permitted emission is centered at 
zero velocity while the forbidden emission is blueshifted along with the wind-indicative He I 1.083 
fim absorption. Variability over the time series is apparent, with the forbidden lines stronger during 
epochs of higher extinction and the H I lines strong at lower extinction; there is little variability in 
the permitted metallic lines. 
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both the optical and near-infrared) in the forbidden line absolute fluxes. Infrared H2 lines behaved 
in a manner similar to the forbidden atomic lines. As the permitted atomic emission is generally 
associated with the accretion flow (with some lines influenced by outflow physics), and the for- 
bidden atomic and the H2 lines are generally associated with outflow, we suggest below that their 
respective correlation and anti-correlation with the photometric state of PTF lOnvg suggests that 
the forbidden and H2 lines are formed outside the extinction region. Other molecular lines, notably 
infrared CO and optical TiO/VO, however, followed the emission behavior of the permitted atomic 
lines, suggesting that these lines are also formed in the inner disk region. 

At high-dispersion, the emission line spectrum of PTF lOnvg bears significant resemblance to 
that of VI331 Cyg in terms of the species present and the relative line strengths. Direct comparison 
of small portions of the red optical and of the 1 jim. region high-dispersion spectra are provided 



in Figure 19 The PTF lOnvg spectrum is more "extreme" in terms of emission line presence and 
strength than the spectra of other well-known "continuum plus emission" sources. Specifically, 
compared to DR Tau, Z CMa, or SVS 13, for example, PTF lOnvg exhibits a larger number of 
emission line species and with larger blueshifts. However, at 1 /jm the blueshifted absorption trough 



in He I A 10830 for PTF lOnvg (Figure 18) is not as deep or as wide as in these objects (see Figure 1 
of Edwards et al.|[2003 for profiles) despite the similar line terminal velocities. V1331 Cyg similarly 
has a He I profile that is much deeper and broader than that of PTF lOnvg. Unlike the spectra of 
the above sources around 1 /um, however, V1331 presents a similarly rich emission line spectrum 
as PTF lOnvg in both the permitted and some of the forbidden line species. In the optical, the 
PTF lOnvg spectrum (from 2010, December, the only epoch in our high-dispersion data set that 
is not limited by signal-to-noise) is again well-matched to V1331 in the permitted line emission 
though PTF lOnvg has stronger and more blueshifted forbidden line emission. 

In the sub-sections below we discuss the details of the molecular, permitted atomic, and for- 
bidden atomic lines and their time series behavior. Summarizing our findings, the evolution of 
PTF lOnvg's spectroscopic properties support our interpretation from the photometry that the 
source undergoes large variations in line-of-sight extinction (A v > 30 m ) towards the origin of the 
optical/near-infrared continuum and permitted emission lines. Because emission equivalent width 
increases (decreases) with continuum fades (rises) as io°- 4x<5m , we do expect some evolution in 
line strength. However, the emission lines which are thought to arise from shocked circumstellar 
material do not appear to be subject to the bulk of these extinction variations; the fluxes of these 
lines are relatively stable throughout the monitoring period, and therefore appear significantly more 
prominent at epochs when the other components of the spectrum appear heavily extincted. 



5.3.1. Molecular Emission 

Above when describing the behavior of PTF lOnvg's optical continuum, we discussed that the 
broad TiO and VO molecular emission at optical wavelengths is time variable. To quantify the 
effect, we measured the strength of the TiO emission at 6250, 6760, 7100, 7800, 8455, and 8880 
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Fig. 19. — Portions of our high-dispersion spectra for PTF lOnvg (in black) compared to the well- 
known emission line dominated object V1331 Cyg (in red). There is a clear spectral similarity 
between these two young stars. Top panel: Keck/HIRES data in the red optical spectral region, 
showing PTF lOnvg in 2010, December. Bottom panel: Keck/NIRSPEC data in the Y-band, 
showing PTF lOnvg on 2010 July 16, near peak outburst. Blue lines identify various permitted 
metallic emission lines, which are similar in the two stars. H I, He I, and forbidden metallic lines 
(located in other spectral regions) are also present in both stars but exhibit different kinematic 
signatures. 
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A as well as the VO emission at 7445 and 7865 using band indices as described in Hillenbrand 



(1997). There is stronger optical molecular emission relative to the continuum during the 2010 



outburst and in early 2012 when the source was bright, relative to the periods when the source was 
fainter, notably in late 2011. This is consistent with the inference drawn from direct examination 



of the spectra both on an absolute scale and on a continuum- normalized scale (see Figure 10). In 
addition to the direct correlation between the strength of the different molecular bands there is 
also a correlation between molecular band strength and continuum slope. 

Molecular emission is also seen in PTF lOnvg's near-infrared spectra from bright photometric 



states, specifically in bands arising from H2O, CO, and VO (Figure 12). Like TiO and VO, these 
bands are commonly seen in absorption in M-type photospheres. Unlike TiO and VO, however, 



CO overtone emission is not uncommon in embedded protostars; for example, Doppmann et al. 
(2005) and Connelley & Greene (2010) detected CO emission from ~15% of the Class I and and 



~25% Class I/II (i.e. flat spectrum) sources in their samples. Detailed analyses of the kinematic 
profiles of these CO bandheads in several particularly well-studied CO emission sources suggest that 
the emission arises from the innermost regions of circumstellar disks in Keplerian rotation (e.g., 
Najita et al.|l996[ ). As Figures [121 and Il7l demonstrate, the strength of the CO and H2O emission 
features are highly correlated with both emission from permitted atomic lines (H I, Na I, Ca I, 
etc.) and with the brightness of the underlying continuum. These molecular band and permitted 
atomic emission lines are all prominent in PTF lOnvg's bright states and all weaken relative to the 
continuum as the continuum fades. 

Summarizing the optical and infrared spectroscopic time series, there are several oxygen species 



molecular lines that are more prominent when PTF lOnvg is brighter. As discussed in Hillenbrand 



et al. (2012) this could be due to molecule formation following dissociation during the bright state, 



or perhaps simple heating and cooling related to variation in optical and higher energy photon flux. 
Alternately, the correlation could be caused by obscuration of the molecular emitting region during 
the faint stages, which is the same explanation we invoke below for the time variable behavior of 
the atomic line emission. 

Molecular emission from H2, however, displays quite the opposite trend: while H2 lines are 
present at all epochs, these lines appear most prominent during PTF lOnvg's faint states. The 
absolute line flux remains relatively constant, however, indicating that the increased prominence 
of these features during faint states is due to improved contrast with the nearby continuum. This 
behavior is qualitatively visible in Figure [l7j where the proximity of the 2.12 /mi H2 line to 
the 2.166 fim H I line provides a clear demonstration that H2 strengthens relative to H I as the 
extinction (as inferred from near- infrared colors) increases. This behavior strongly suggests that 
PTF lOnvg's variable extinction obscures the emission region responsible for the H I and other 
permitted emission lines, but not the emission region from whence the H2 emission arises. To test 



this expectation in more detail, we compare in Figure 20 the ratio of the H I and H2 line strengths 
to a simple model in which H2 line strength is invariant while H I line strength diminishes. This 
simple model accurately replicates the dependence of the H2-to-HI line ratio on Ak for all but the 
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Fig. 20. — HI 2.16 /im - to - H2 2.12 fj>m line ratio measured from moderate resolution near-infrared 
spectra (asterisks) as a function of PTF lOnvg's K s magnitude relative to peak brightness. Error 
bars for the ratio measurements were established by adding synthetic noise to the spectrum then 
re-measuring the line ratios. Overplotted for comparison as a dashed line are the predictions of 
a simple model in which the H2 line flux is constant with time, while the HI 2.16 ^m line flux 
diminishes as expected if the differential K s magnitude is due entirely to extinction (i.e, HI = 
HI c x io~°- 4xA ^ s ). The anomalous point at high extinction corresponds to the 2011, July 15 epoch 
where the brackett gamma line strength is very weak, and therefore may be more susceptible to 
systematic errors introduced by mis-cancellation of the HI profile of the telluric standard. Except 
for this anomalous point, the observed line ratios agree quite well with the quantitative predictions 
of the extinction variation model. 
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most heavily extincted epoch (for which increasing systematic errors may lead to underestimated 
errors based on photon statistics alone). 

Additional insight into the physical processes responsible for the H2 emission can be extracted 
from the moderate resolution near-infrared spectrum obtained on 2011 July 15, when PTF lOnvg 
was near a photometric minimum. This spectrum is nearly featureless, save for a number of well 
detected H2 and [Fe II] emission features. This provides a unique opportunity to characterize the 
physical conditions responsible for the H2 emission. Shocks due to stellar winds and outflows, 
as well as photo-excitation by high-energy radiation, are both capable of producing H2 emission. 
Previous studies of Classical T Tauri stars and Class I protostars have identified systems that 
appear to be templates for each type of H2 excitation mechanism, as well as numerous systems 
where the excitation mechanism remains ambiguous ( Beck et al.||2008 Greene et al. 2010). 



Considering first the 2010 July 14 spectrum, Covey et al. (2011) were able to confidently 



measure only the 2.12 /jm H2 line, as well as state a tentative detection of the 2.24 ^m H2 line. 



From the ratio of these two lines, Covey et al. concluded that shock excitation was the most 
likely mechanism for exciting PTF lOnvg's H2 emission, but noted that this conclusion remained 
uncertain. The 2011, July 15 spectrum, by contrast, enables clean measurements of nine K-band 
H2 emission lines, whose strengths we present in Table |8j For comparision we also provide the line 
strengths expected for different emission mechanisms: first as predicted by models of J- and C-type 



shocks in stellar outflows (Smith 1995) and by models of UV and/or X-ray photo-excitation in 



the presence of circumstellar dust (Nomura et al. 2007), and second, empiricially measured values 



for spatially extended H2 emission near RW Aur (Beck et al. 2008) and for a knot in the HH 54 



outflow ( Giannini et al.||2006 ). Finally, Table [8] also summarizes the observed and predicted values 
for specific line ratios identified by Beck et ah] ( |2008 ) as particularly useful for disentangling the 
underlying excitation mechanism for H2 emission. The PTF lOnvg line ratios are extremely similar 



to those measured by Beck et al. (2008) using adaptive optics imaging that spatially resolves the 



circumstellar H2 emission near RW Aur. The line strengths measured by Giannini et al. (2006) 



from Knot B of HH54 are also in reasonable agreement with our observations of PTF lOnvg. This 
similarity is consistent with the hypothesis that the origin of PTF lOnvg's H2 emission is a region 
that is well separated from the central protostellar source. Models that reproduce the line ratios 
observed for PTF lOnvg and RW Aur include those calculated by Smith (1995) for a ~2000 K 



C-type shock and the models produced by Nomura et al. (2007) for pure X-ray photo-excitation. 



A final clue to the origin of PTF lOnvg's H2 emission is provided by the kinematics of the 
lines, which are uniformly blueshifted by ~200 km/s with respect to the rest velocity of the North 
America Nebula as well as the observed velocities of PTF lOnvg's permitted atomic features. 
This kinematic profile strongly suggests that PTF lOnvg's H2 emission originates from within a 
stellar jet; interestingly, isolating the location of this emission region does not completely resolve 
the ambiguity between X-ray and shock heating as potential excitation mechanisms, as Chandra 
observations of the HH 154 protostellar jet revealed that shock regions can generate X-rays in situ 
in emission regions with expansion velocities of ~ 500 km/s. Figure 17 demonstrates the behavior 
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of H2 (second panel from right) relative to [Fe II] (third and fifth panels from right) which vary 
in a correlated manner over the time series, and relative to CO (right-most panel) and various 
permitted atomic lines (all panels) where there is anti-correlated behavior. 

In summary, multiple lines of evidence suggest that PTF lOnvg's H2 emission arises from a 
spatially extended outflow. First we note the insensitivity of the emission strength to the variable 
extinction that apparently occults the innermost disk and accretion region where the near-infrared 
continuum and the other molecular (see above) emission as well as the permitted atomic (see 
below) emission lines form. Second we note the similarity in H2 line strengths to RW Aur's spatially 
extended emission and to models of excitation due to mechanical shocks and X-ray photoexcitation, 
both of which have been associated with other protostellar jets. 



5.3.2. Permitted Atomic Emission 

Over most of the time series, the strongest emission lines in the optical spectra of PTF lOnvg 
are those of Ha and Ca II, a typical situation for young stars. As discussed below, during fainter 
photometric stages, however, the forbidden lines became stronger than these atomic lines. The Ha 
and Ca II line strengths vary both absolutely and relative to the red optical continuum (top and 



bottom panels of Figure 10). 



The Ca II triplet emission varied between —15 and —25 A but became much weaker in 2011, 
June, when the continuum faded. For classical T Tauri stars in Taurus, the line strengths among 
the three lines of the Ca II triplet are typically 1.20 ± 0.19 for the 8542 A to 8498 A ratio and 
1.28 ±0.13 for the 8542 A to 8662 A ratio (based on analysis of our own collection of optical spectra 
for representative samples). In PTF lOnvg, however, the 8498 A line is particularly strong; the 8542 
A to 8498 A ratio has a mean and dispersion of 0.93 db 0.15 over the time series (considering only 
spectra for which the lines are measureable), while the 8542 A to 8662 A ratio is consistent with 
expectations at 1.14 ± 0.18. Although highly unusual, these observed line ratios for PTF lOnvg are 



not unique. Hamann & Persson (1989) found for V645 Cyg ratios among the triplet lines of 0.85 
for 8542 A to 8498 A and 1.1 for 8542 A to 8662 A, similar to those in PTF lOnvg. Notably the 
line ratios are far from their expected ratio of 1:9:5 for optically thin emission, which is typically 



the case for young accreting stars. Kwan Sz Fischer (2011) find that the physical conditions for 
producing strong Ca II triplet emission include relatively high density njj ~ 10 12 cm~ 3 and low 
temperature T <7,500 K, with the line formation taking place in the accretion flow. 



As illustrated in the high-dispersion profiles of Figures 15 , the Ca II triplet lines all peak to 
the red of zero velocity and the velocity profiles are similar between the 2010, December and 2011, 
December spectra (when atomic emission was prominent). However, there are also indications of 
differences among the three lines of the triplet. The 8498 A line is the broadest of the three, 
with FWHM~200 km/s, has an approximately symmetric profile, and exhibits no evidence for 
sub-continuum absorption. The 8542 A and 8662 A lines, by contrast, are somewhat narrower 
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and asymmetric (narrower on the blue side than the red side), with sub-continuum "dips" in their 
blueshifted components that change between the two high dispersion spectra and indicate wind 
absorption. The 8542 and 8662 A structure is similar to that seen in other strong outflow sources 
such as V1515 Cyg and V1057 Cyg. On its red side, the 8498 A line traces well the 8662 A profile 
but the 8542 A line is broader, on the blue side, though, it is the 8542 and 8662 A profiles that are 
quite similar. While these three lines share the 2P— >2T) transition, they differ in their spin with 
the 8542 and 8662 A electrons changing spin between the 4p and 3d energy levels, and the 8498 A 
electrons retaining the same spin, which is perhaps a clue to the apparent differences in line optical 
depth. 

In Ha, the emission line equivalent widths varied at low dispersion between —20 and —60 A 
during brighter photometric states but became quite large (many hundreds of A) as the source 
continuum level faded during late June through August of 2011. The strength and profile behavior 



as discerned from our collection of low resolution spectra is illustrated in Figure 13 while the 



high-dispersion behavior is illustrated in Figure 16 



The low resolution spectra from 2010 July, August, and September all show a redshifted 
emission shoulder in Ha, with the first July spectrum exhibiting a clear blueshifted absorption 
component, which is then ambiguous in the August and September data. The line profile as seen 
in the high-dispersion 2010, December spectrum shows more clearly both the redshifted emission 
component and a two-component structure to the blueshifted absorption, similar to the structure 
seen in the Ca II lines described abov^J In H/3, the line structure has the same —250 km/s 
absorption component as seen in the Ha line, but a deeper —125 km/s absorption component that 
extends to lower velocities and even to the positive side of the line. During the faint photometric 
state spectra of 2011, June and 2012, January, the Ha profile lacks any redshifted emission and only 



weak blueshifted emission is observed (see Figure 16 ) while the brighter state 2011, December profile 



is similar to the 2010, December (bright state) profile with similarly strong blueshifted emission 
though much weaker redshifted emission. The blueshifted emission components in H/3 and Ha are 
similar and exhibit similar time series behavior. Specifically, while from 2011, May through to 



2011, October, the Ha line is considerably weaker (Figure 10), the even greater dimming of the red 
continuum produced enhanced line equivalent width during this period (left panel of Figure 13). 
The 2011, October and November profiles then revert back to the redshifted emission shoulder, as 
the source brightened. 

This time series behavior suggests that at least some of the blueshifted Ha emission originates 
in the outflow or outflow shock region, similar to the H2 emission (dicussed above) and to the 
optical/infrared forbidden line emission (discussed below), both of which are similarly blueshifted. 
The blueshifted velocities for all of these emission tracers at the position of the continuum source 



Based on optical spectra taken on 2010 September 5 and 2010, November 25, 



Aspin 



(2011) claim that the 



blueshifted Ha absorption reported by Covey et al. (20111 had disappeared. Our data show that the absorption 
weakened in the low dispersion data but it clearly remained present. 
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are similar to the velocities reported in section 5.3.4 for the spatially offset blueshifted emission 
that is observed, and arises from the interaction of the outflow with either ambient material in the 
local vicinity or itself. 

In both the Ca II and the H I profiles there is no apparent change over the existing time series 
in the velocity of the emission line peaks (which, as described above, are slightly to the red in both 
of these lines) or in the terminal velocity of the blueshifted absorption components. That there is 
little change in the terminal velocity of the wind means that the launching velocity is not changing, 
even though there is changing structure in the absorption part of the P Cygni profiles. 

Notably, no optical He I A 5876 emission is seerj^J even though this line is near-ubiquitously 
detected for classical T Tauri stars in high-dispersion spectra. For strong emission at 5876 A, 
T«I5,000 - 20,000K and n H > W u would be required dKwan & Fischer 112011 1) and thus a lack of 



5876 A emission while lower energy 10,830 A emission is observed could indicate an upper limit on 
temperature and/or density in the emission region. 

Other prominent atomic emission lines in the optical include those of O I, Fe I, Fe II, Ca I. 



Those lines that are visible in the existing high-dispersion spectra (see Figure 15) are generally 



symmetric and centered at zero velocity, as well as narrow in velocity with full- width half-maximum 
~45-50 km/s (e.g. Fe I lines) to ~65-70 km/s (e.g. Fe II lines). Their origin may be on the stellar 
surface or in the disk. The O I 8446 A line, however, though also centered at zero velocity, is 
broadened with full-width half-maximum ~250 km/s. This width is similar to that exhibited by 
the Ca II and Ha lines, and suggests origin in the inflowing and/or outflowing material. In the low 
dispersion spectra we measure a ratio of Ca II 8498 A to O I 8446 A (for which the time variable 
extinction effects are not a factor given the small wavelength difference) that does change somewhat 



over time but notably always exceeds 10. Kwan & Fischer (2011) find that values larger than 7 
for this ratio indicate temperatures T <7500 K and densities nu > 10 12 cm -3 for the emitting gas. 
Only the most extreme accretors among well-studied Taurus stars have similar Ca II 8498 A : O I 
8446 A line ratios. 

Another notable atomic line seen in emission is Li I A 6707. Interpretation of the profile (Figure 
15) should take into account that there is a y 3 F°-e 5 P term transition of Fe I reported by the NIST 



Atomic Spectra Database Ralchenko et al. (2011) at similar wavelength, specifically 6707.43 A. No 



other lines of this Fe I multiplet are within our spectral range and so the possible identification as 
such can not be supported. We note, however, that another Li I line at 6103 A may also appears 
weakly in emission in our spectra, possibly supporting the Li I emission scenario. Although this 
line also has a nearby Fe I potential identification, the w 5 T)° -5/2[7/2] term transition at 6103.64 
A, this is a much higher excitation Fe I line than those found otherwise in the PTF lOnvg spectra. 



Redward, as illustrated in Figure 19 the atomic line emission apparent in the high-dispersion 



A narrow 6678 A emission line might be associated with He I but the absence of the related 5876 A line suggests 
that this is more likely Fe I A 6678. 
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1 /um spectra includes zero velocity C I, S I, Si I, Sr II, and Ti I multiplets. There is no Fe I or Fe II 
apparent. The three epochs of NIRSPEC spectra, which cover dramatically different brightness 
stages of PTF lOnvg, show little difference in the metallic line emission profiles (e.g. Figure 18 ) 
which generally match quite well over a range of line excitation potentials. Changes over time are 
apparent, however, in the H I and He I profiles, with both lines exhibiting blueshifted absorption 
components as discussed in a previous section. Other prominent near-infrared emission lines seen 
in the low dispersion spectra include Na I, Mg I, and Fe II. As Figure 17 demonstrates, these lines 
appear to vary in concert with one another, and are anti-correlated with the extinction estimates 
inferred from PTF lOnvg's near-infrared photometry. 

A common practice is to convert infrared hydrogen line fluxes into mass accretion rates on 
to the central star. Following Muzerolle et al. (1998) and assuming a mass of 0.5 M and radius 
of 3 R.0 for the star, we derive accretion rates of 3 x 10~ 7 MQyr _1 in the brighter states and 
6 x lO 8 M0yr _1 in the fainter states using Br7 with an upper limit of 5 x lO _12 M0yr _1 when 
the line was undetected during summer of 2011. Using Pa/? the values are about a factor of 4 
lower. No extinction corrections are applied in these calculations, and thus the numbers are clearly 
lower limits in the faint state and probably the bright state as well. In addition to the extinction 
uncertainty, it is not clear that the accretion model on which the line luminosity to accretion rate 
calibration is based applies to PTF lOnvg. 



5.3.3. Forbidden Atomic Emission 



In young stars, forbidden lines of e.g. [O I], [N II], [S II], and [Fe II] in the optical and 
near-infrared spectra are attributed to few hundred km/s radiative shocks (those that cool by such 
radiation more quickly than they would by expansion) in outflowing material in jets associated with 
accretion disks. Unlike the permitted optical/infrared emission lines discussed above, the forbidden 
emission lines are (almost always) optically thin. They trace plasma over a wide range in density 
n e ~ 10 3 — 10 8 cm -3 in the temperature range T ~5, 000-20, 000 K - for hydrogen ionization fraction 
n e /n}{ ~ 1/10, which is appropriate for the mostly neutral young star jets. 

The moderate temperature, low density forbidden lines are formed closest to the shock, with 
temperature decreasing and density increasing further from the shock in the pre-shock medium. 
The infrared H2 lines discussed above are generally thought to form further away from the shock 
front than the optical/infrared forbidden lines, probing cooler temperatures in either post-shock 
gas along the jet axis that has sufficiently cooled, or material perpendicular to the jet axis that is 
entrained/stirred by the flow. As with the density and temperature gradients, the ionization fraction 
is likely not constant along the flow. The implied range of physical parameters in the outflow 
sampled in spectroscopic observations such as ours that do not spatially resolve (see, however, 
section 5.3.4) the cooling region/s in the outflow zone means that we should consider the forbidden 



line emission indicative but not diagnostic of jet/outflow conditions. See Hartigan et al. (2000) for 



a review and Brugel et al. (1981) for a classic study of forbidden lines in Herbig Haro objects. 
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Indeed, the typical [O I], [N II], [S II], and [Fe II] forbidden lines are present in PTF lOnvg spec- 



tra, but display significant variability over time (see Figure 16). The forbidden profiles are purely 
blueshifted, with a range of profile shapes. Notable in the spectroscopic time series is the increasing 
prominence of these and other forbidden lines in 2011, June - as the broadband optical/infrared 
continuum faded to below detectable levels in our PTF R-band photometry. Specifically in the red 
optical spectra, while the somewhat common forbidden line doublets of [O I] 6300, 6363 A, and 
[S II] 6717, 6731 A as well as [Fe II] 7155, 7172 A and [Ca II] 7291, 7324 A, had been weakly present 
in the 2010 spectra and in the first 2011 (May) spectrum obtained, by 2011, June these lines had 
signficantly increased in strength relative to the continuum. Over just a three week time frame, 
they along with other somewhat unusual forbidden lines (such as the [C I] 9824, 9850 A doublet, 



[Cr II] and [Ni II]; see Figures 10 and 13 14) rose in prominence to dominate the optical spectrum. 



Notably, emission in [N II] is seen only during the very faintest epochs. This increasing prominence 
of the forbidden lines was accompanied by a weakening of the permitted emission lines, most visibly 
the Ca II triplet and several Fe II lines that previously were quite strong. This behavior suggests 
that we are seeing, in the faint photometric states, emission from a Herbig-Haro flow that is not 
visible when the source is brighter. 

We can gain insight into the physical conditions necessary for these forbidden line observations 



by taking guidance from the material assembled in Dougados et al. (2010) and references therein. 



For the different forbidden line species we discuss first the line strengths and their variation, and 
then the line profiles and their variation. 

For the [O I] 6300 and 6363 A doublet, having the highest critical density of the red optical 
[O I], [N II], and [S II] lines, the evolution of the observed line ratio is truly remarkable. Typically 
observed [O I] 6300 A : [O I] 6363 A values are in the range 2.5-4 : 1 with the theoretical line 
ratio 3:1 based on Einstein A values. Before and after the deep photometric fading episode of 2011, 
the expected 3:1 ratio for optically thin emission was indeed consistently observed to within 10% 
in PTF lOnvg spectra. However, as the continuum flux faded and the forbidden line prominence 
increased, the [O I] doublet ratio was variable between the typical 3:1 and values as extreme as 1:2 



in several of the faint state 2011 June, July, and August low dispersion (see Figure 13), as well as 



high-dispersion (see Figure 16) spectra. The variability in the doublet line ratio suggests dramatic 



changes in the optical depth of the [O I] lines, from optically thin to a likely optically thick case 
where 1:1 is expected; we can not, however, explain the observed situation that the 6363 A line is 



stronger than the 6300 A line at some epochs. We note that Li & McCray ( 1992 ) developed a model 
for [O I] 6300 A : [O I] 6363 A doublet emission in SN 1987A, which exhibited a systematically 
time variable ratio that increased from an optically thick 1:1 to the usual optically thin 3:1. The 
model was based on an expanding shell with a small filling factor, i.e. clumpy distribution of 
[O I], that is not necessarily applicable here given differences in expected geometry between a 
spherically expanding supernova shell and a more collimated young star outflow, but does suggest 
that variable line optical depth is a plausible explanation for our observations. Notably, the [O I] 
changes observed in PTF lOnvg were not gradual and systematic as observed in SN 1987A; rather, 
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they were abrubt with the 6300 A : 6363 A ratios changing literally from night to night, e.g. from 
1:1 on 1 August, 2011 to 3:1 on 2 August, 2011 when the data derive from the same instrument, 
configuration, observers, slit position angle, etc. so the variability can not be attributed to spatial 
alignment or other similar obserational effects. Plausible explanations for our observations, whether 
related to optical depth effects or not, need to account not only for the truly peculiar 1:2 line ratio 
but also the short term variability of the ratio. 

At all epochs, including during the photometric fading period when the [O I] ratio is out of the 
expected 3:1 ratio, the much lower critical density [S II] 6717, 6731 A doublet notably remained 
within the expected range for its ratio. Typically observed values in young stars are 0.5-1 : 1, with 
the theoretical ratio in the high density limit 0.43. The PTF data are consistent with these values 
but with some variation over the time series. Most of the data exhibited line ratios close to the 
high density expectation (indicating n e > 10 4 cm -3 ) but some epochs show more equal strength 
doublet lines, indicating lower density gas (n e ~ 10 3 cm -3 ). 

Notably, neither [O III] 4959 / 5007 with a similar critical density to the strong [O I] doublet, 
nor [O II] 3726 / 3728 with a similar critical density to the strong [SII] doublet, are seen. Although 
our spectra become somewhat signal-to-noise challenged in the blue, strong line emission compara- 
ble to the red line emission would be detectable. The [O III] line absence provides a constraint on 
the electron temperature of the emitting gas and is, indeed, rarely seen in well-studied HH objects. 
The [O II] absence can be taken to confirm the assumption of a low ionization fraction. 

The [N II] 6548, 6583 A doublet is of intermediate critical density and has a typically observed 
ratio of 1-5 : 1 in young stars with a theoretical ratio of 3 based on Einstein A values. These lines 
are present in just a few of our low resolution spectra during the faintest photometric states, but 
are hard to disentangle from the broad Halpha emission (see Figure 13). They were only weakly 
visible in our high-dispersion data. We note that [N II] 6548, 6583 A emission is seen, however, at 



positions that are spatially offset from the PTF lOnvg point source, as discussed in section 5.3.4 



Some of the weak emission seen in Figure 16 may arise in the extended emission region closest to 
the trace (see Figure 21). 

Considering a higher density diagnostic, the infrared [Fe II] 1.533 and 1.644 (im lines, like the 
optical forbidden lines, are more prominent when the source is faint in the continuum. The ratio 
of these two lines is density sensitive and indicates n e > 10 5 cm -3 . The 1.257 and 1.644 /im line 



ratio can be used to measure extinction since the lines share the same upper level (see Connelley 



Sz Greene (2010)). For PTF lOnvg the time series data indicate low extinction (Ay < 2 mag) to 



the jet region - in contrast to our earlier findings on the large and time variable extinction towards 
to near-infrared continuum. The above and other infrared Fe II lines that are prominent during 
PTF lOnvg's faintest epochs are tabulated in Table [8] and can be compared in line strength ratio to 
that observed by Giannini et al. (2006) from Knot B of HH54. The two sets of measurements agree 
moderately well, providing additional support for the existence of a significant outflow within the 
PTF lOnvg system. The optical [Fe II] lines probe even higher density material. The measured 7155 
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A, 7452 A, 8617 A, and 8892 A lines of [Fe II] yield ratios that according to Bautista & Pradhan 
(1998) indicate gas at densities of n e ~ 1 — 5 x 10 6 cm -3 in the vicinity of PTF lOnvg. The optical 



[Fe II] doublet at 7155, 7172 A varies in our time series, strengthening during the 2011 faint epochs 
but perplexingly failing to re-emerge in the late 2011 / early 2012 fade like other forbidden lines. 
These lines also exhibit a ratio both well above and below its expected 3:1 ratio for optically thin 
emission. Unlike the time series behavior of the [O I] doublet ratio, however, the [Fe II] anomalies 
are not restricted to the faint epochs. There is no correlation of the [Fe II] line ratio variations and 
the [O I] line ratio peculiarities. 

High density gas is also probed by [Ni II] and [Ca II] lines. The observed ratio of [Ni II] 7412 



A to 7378 A in the 0.3-0.5 range indicates densities n e > 5 x 10 cm 3 according to |Bautista et 



al. ( 1996 ). The observed ratio of Ca II 8542 A to [Ca II] 7291 A varies between 8 and 15 during 
the bright photometric states and is indicative of densities in the range n e ~ 3 — 5 x 10 8 cm . 
During the faint states this ratio becomes very small as the atomic emission is not detected while 
the forbidden emission becomes quite strong relative to the continuum. 



As illustrated in Figure 19, there is also forbidden line emission in the high-dispersion 1 /im 
spectrum, with a blueshifted set of [N I] doublets around 1.040 /tm (specifically 1.040059, 1.040100 
and 1.041002, 1.041044 /mi) and a blueshifted [S II] series located at 1.029, 1.032, 1.034, 1.037 
/im. They are indicative of hotter temperature gas, in the 30,000-40,000 K range, than the optical 
forbidden lines which probe 5,000 to 20,000 K temperatures. Like the optical lines, the near-infrared 
forbidden line strength appears to correlate with the photometrically derived extinction estimates 
(Figure [9]); specifically, the 2012 spectrum exhibits stronger [N I] and [S II] lines than the 2011 
and 2010 spectra, both of which were taken when the extinction was about half of the value for 
the 2012 spectrum. The non-simultaneous optical-infrared spectroscopy in combination with the 
unknown extinction and evidence for large extinction variations over time precludes us from using 
infrared-optical line pairs of these species to probe either temperature (e.g. from optical 6730 A vs 
infrared 1.03 /jm [S II] at high temperature, or optical 8620 A vs infrared 1.64 /mi [Fe II] at lower 
temperature) or ionization fraction (e.g. from optical [N II] and infrared [N I]). We simply note the 
presence of these hotter lines in the 1 /im region spectra. 

In addition to the line strength variations, there were variations in the line profiles of the 
forbidden emission in PTF lOnvg. These are illustrated by our high-dispersion optical data as 



shown in Figure 16 In all forbidden lines and at all epochs, the profiles are purely blueshifted 
with maximum emission velocities of ~250 km/s. At the earliest 2010, December epoch, the [O I] 
doublet has a broad FWHM and nearly square profile ranging from to —250 km/s, peaking at 
about —130 km/s, whereas at the later epochs the profile becomes skewed with a broad base but 
increasing flux towards higher velocities, peaking near the terminal velocity which remains at the 
same —250 km/s terminus. The skewed profile is present both in fainter and brighter continuum 
states of PTF lOnvg. Similar line profile morphology and behavior over time is seen in the [Fe II] 
doublet. The [S II] doublet seems to have the opposite trend, however, with a similarly shaped 
but narrower and slightly lower velocity (—100 km/s) peak relative to the [O I] doublet early on, 
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but becoming more boxy and broad as the [O I] became more skewed. The [N II] doublet appears 
quite similar to the [S II] doublet but is much weaker, and essentially undetected in most of our 
high-dispersion spectra. 

The [Ca II] doublet, unlike the above lines, is double-peaked, with a higher velocity component 
centered at about -150 km/s and extending to —200 km/s (similar to the [S II] and [N II] doublets), 
as well as a second, nearly zero velocity but slightly blueshifted peak. This double-peaked profile 



(see Figure 16). is seen in the 2010, December and 2011, December spectra, both taken near 
photometric peaks in the lightcurve. During the fainter photometric states, however, only the 
blueshifted peak is seen, and presenting the same skew presented by the [O I] and [Fe II] doublets. 
A similar double-peaked and time-series behavior is exhibited by the optical [C I] lines. The [Ca II] 
doublet, unlike other forbidden lines, has a possible origin in photon pumping through the permitted 
Ca II 8498, 8542, and 8662 lines rather than requiring collisional pumping (though collisions may 
still be the dominant mechanism in young star jets). The 7324 A: 7291 A line ratio of [Ca II] is 
consistent with the expected equilibrium value of 0.67. 

In summary, the forbidden emission line strengths in PTF lOnvg indicate a wide range of gas 
densities, from n e ~ 10 3 cm -3 to n e > 10 8 cm -3 . Line strength variations some line ratio variations, 
and line profile variations were all seen over the time series. 



5.3.4- Spatially Offset Emission 

The velocity structure along the spectral trace as defined by the (sometimes weak) spectral 
continuum exhibits blueshifted emission in a large number of forbidden lines as well as in Ha, 
notably during the fainter photometric states of PTF lOnvg's lightcurve. This on-source emission 
was discussed in depth above. Many of the HIRES images additionally exhibit blueshifted emission 
that is spatially offset from the continuum source position. This is illustrated via position-velocity 
diagrams in the vicinity of Ha and forbidden [O I] 6300 A, [N II] AA 6548, 6583, and [S II] AA 6716, 



6730 in Figure 21 



We show the two-dimensional HIRES spectra at three epochs: 2010 December when PTF lOnvg 
was relatively bright, 2011 June when it was quite faint, and 2011 December when it was again 
bright. The two December spectra were obtained at similar position angle (PA), 110.8° and 102.6° 
for the 2010 and 2011 observations, respectively, with a 7" long decker, while the 2011 June spectrum 
was obtained with a 14" long decker and a PA of 176.0°. Despite the differing PAs between the 
spectra, there are a few features that are commonly observed at all epochs. These include a bright 
knot of emission extending ~2-6" from the stellar trace at ~ —190 km s , which can be seen in all 
epochs in Ha, [S II], and [N II] A 6583 A (note that only the 2011 June spectrum had a long enough 
decker to see the emission 6" from PTF lOnvg). This same knot is also visible in [N II] A 6548 A in 
the two 2011 spectra as well as [O I] A6300 A in the 2011 December spectrum. In addition, there 
is plausible but weak offset emission in [Ca II] 7291 A at this same —190 km/s velocity, but in the 
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Fig. 21. — Spatial-spectral diagrams centered around the optical Ha, [O I], [N II], and [S II] lines. 
Units on the ordinate are arc seconds while those on the abscissa are given in both wavelength (A) 
and velocity (km/s) with each panel showing a 350 km s _1 portion of the emitting region. The 
image intensity scaling is linear with the range set over a range of counts selected for each feature 
to emphasize the spatially extended emission (not the same in all panels). Cosmic rays have not 
been removed from the images. "Sky" emission fills the decker vertically at low velocity locations 
that depend on heliocentric correction and may not be the same from line to line depending on 
whether the emission is dominated by telluric vs nebular background emission. The left and right 
panels show the HIRES observations from 2010, December (PA = 110°. 8) and 2011, December 
(PA = 102°. 6), respectively. The middle panel shows data from 2011, July (PA = 176°. 0) when 
the decker setting allowed twice the spatial coverage of the other two panels but otherwise exhibits 
similar structure away from the main trace. Blueshifted emission is seen below the trace at several 
velocities (roughly -190, -160, -140, and -20 km/s in all columns plus at -75 and blueward in the 
middle panel only, with its larger extent in the spatial direction). Blueshifted emission is also seen 
above the trace at about -50 km/s. 
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7324 A line of the doublet it is not obviously offset from a sky line that coincident ally appears at a 
similar velocity; these images are not shown. Fainter knots, apparently extending from the position 
of the continuum source, can be seen a few arcsec below the trace at ~ —140 and —160 km s _1 . 
These features are present in all three epochs and are most easily seen in [S II]. The final knot that 
can be seen in all epochs is ~2" above the stellar trace at ~ —50 km s -1 , and it too is most readily 
identifiable in [S II]. 

Interestingly, in the 2011 June spectrum there is virtually no [N II] emission from the PTF lOnvg 
continuum position, and for A6583 in particular, it can be seen that each of the aforementioned 
knots are connected spatially. This structure is clear evidence that PTF lOnvg is driving a jet with 
a terminal shock ending ~4" from the star and maximum velocity ~200 km s _1 . The appearance of 
these features at all epochs, regardless of the PA, suggests that PTF lOnvg is driving a wide-angle 
outflow. 

There are additional spatially offset features which are not observed in each epoch. The most 
prominent is a knot seen in the 2011 June spectrum at ~ —20 km s _1 located ~4" below the trace 
of PTF lOnvg and extending to the edge of the decker ~7" from PTF lOnvg. This feature is most 
easily identified in [O I] since it is blended with diffuse nebular emission in all the other panels. This 
feature is too far from PTF lOnvg to be detected in the 2010 December spectrum as the decker did 
not extend 4" below the PTF lOnvg continuum, while the same emission feature may be present in 
the 2011 December spectrum, though blending with nebular emission (or the [O I] night sky line) 
make a positive identification challenging. Finally, the last prominent feature, which is seen only 
in the 2011 June observation, is a point source-like trace seen ~7" below PTF lOnvg. This trace 
is seen m Ha, [S II], and [N II] A 6583 A. Optical imaging does not reveal any point sources within 
~8" of PTF lOnvg and so the origin of this emission is currently unclear. High-angular resolution 
narrow-band imaging of the field around PTF lOnvg should elucidate if this feature is from another 
star or somehow related to PTF lOnvg itself. 

No spatially offset emission is apparent in the NIRSPEC 1 /im high-dispersion spectral imaging 
data, e.g. in [N I] or [S II], or in the 1.083 /im He I line. For these observations the slit length 
was 12" and thus able to sample the relevant spatial range where offset emission is seen in optical 
forbidden lines, as discussed above. 

6. Constraints on Multiplicity 

Neither the raw nor the reduced and coadded adaptive optics imaging frames from Keck/NIRC2 
show evidence for point source companions within several arcseconds of PTF lOnvg. Figure [22| 
shows the sensitivity of the observations to off-axis sources. The diffraction-limited data rule out 
the presence of companions down to a level of AK S = 3.2, 5.0, 7.0, 7.4 mag at 0.25" , 0.50" , 1.0" , and 
2.0" (5a) respectively. For comparison, the typical FWHM of our seeing-limited PTF observations 
in the R-band and the beam size of our CARMA data at 2.7 mm are both 2.0". 
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Fig. 22. — Contrast versus angular separation curve for the field centered on PTF lOnvg. A 
maximum contast of 7.5 mag is achieved in the K s band beyond 1".0 separation, corresponding to 
520-600 AU depending on the adopted source distance. Detailed image analysis reveals that no 
stellar companion candidates are detected above these limits though the observations can not rule 
out very low mass stellar or any substellar companions. 
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With a primary source brightness of K~8.5 mag at the time of the AO observations, the 
achieved contrast limits correspond to a lack of companions brighter than K~16 mag outside of 
the optical and mm PSF. At closer separations, the flux limits can be translated into upper mass 
limits to any compansions of 1.14, 0.66, and 0.30 Mq at separations of 130, 260, and 520 AU. These 



numbers use the isochrones of Girardi et al. (2002) and assume an age of 1 Myr, a distance of 520 



pc, but do not include a foreground extinction estimate since this is unknown (but would raise 
the mass limits). The companion upper mass limit sensitivity is not particularly deep despite the 
large contrast achieved with Keck/NIRC2 since the 2 /xm brightness of PTF lOnvg is dominated 
by circumstellar excess rather than pure photosphere. 

For reference, companions to several well-known FU Ori-type stars do reside within this sep- 
aration and contrast range range: at 196 AU and AK S = 4.1 mag for FU Ori (Wang et al. 2004 



Pueyo et al.]|2"012"] ) and at 115 AU and AK S = 1.2 mag for Z CMa ( |Koresko et"aTfl99l) |Hinkley et 



al. 2012). Although the companions have been claimed to have some relevance to the outbursting 



behavior of FU Ori stars, the ability for such wide companions to instigate the inner disk instability 
is limited. Companions are also invoked in several examples of presumed disk-eclipsing young stars 



that, like PTF lOnvg, have long periods. Specifically, KH15D at 48 days (Hamilton et al. 2012), 



WL4 at 131 days (Plavchan et al. 2008), and YLW16A at 92 days (Plavchan et al. 2010) have 



extinction-like flux dips between a few tenths and several magnitudes in depth. Several of these 
stars have companions in the 25-50 AU region that are hypothesized to induce warps in the disks. 
Our existing AO observations are not sensitive to such close binaries. 



7. Discussion and Interpretation 

7.1. The Nature of the Large- Amplitude Variable PTF lOnvg 

In this section we provide broader context for the observations of PTF lOnvg that we have 
described above, namely: 1) the large amplitude optical/near- infrared magnitude and color vari- 
ations that are broadly consistent with expectations from extinction variations; 2) the periodic 
nature of the multi-band lightcurves; 3) the respective correlation / anti-correlation of permitted / 
forbidden atomic spectral line emission strength with continuum brightness; 4) the rare detection 
of molecular TiO/VO emission at optical wavelengths, with corresponding molecular CO, H2O and 
VO emission in the near-infrared; 5) the exhibition of strong P Cygni line profiles in several atomic 
species that are sensitive to winds and outflows; 6) the lack of evidence for stellar multiplicity; 7) 
the Class I-type spectral energy distribution; 8) the spatially extended molecular outflow. 



Before the 2010 brightening that was detected by both Itagaki et al. and |Covey et al. the 



source discussed here as PTF lOnvg was known as IRAS 20496+4354 but had not been well-studied. 



Although there was no previously reported near-infrared or optical photometry for the object, Aspin 



(2011) and we in this paper have demonstrated through investigation of archival images that the 



source has been sporadically visible in several past observations. The recent photometric behavior 
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included brightening on a several month time scale followed by an overall fade that has been 
punctuated by semi-periodic spikes in brightness over the past several years; each of the brightness 
peaks subsequent to the absolute peak in 2010 has been fainter, however. 

Both previous to and subsequent to its 2010 brightening, the source can be characterized as a 
Class I (disk-plus-envelope) young stellar object. Spectroscopically, PTF lOnvg is highly unusual 
for a Class I, or any other Class, young star. In addition to the rare TiO/VO emission and more 
common CO emission as well as permitted atomic emisison during bright photometric states, in 
its faint photometric state the source is dominated by forbidden emission lines, some of which are 
commonly seen in young stars but others of which are rarely seen. 

One interpretation is that the object could be in its initial envelope-clearing stage such that the 
embedded optical/infrared point source is in the process of becoming revealed. The consistency of 
both the photometric and spectroscopic behavior with variations in line-of-sight extinction supports 
this interpretation. Further, the molecular emission could be the result of heating and cooling at 
high latitudes in the outflow cavity as stellar photons first reach these areas of the circumstellar 
environment. As an analogy, classical novae, which are also driven by accretion (on to a white dwarf 
of material from an unseen companion) frequently show "dust dips" in their outbursting light curves 
with several episodes of fading and re-brightening ocurring on tens to hundreds of days time scales, 
similar to the situation we have found for PTF lOnvg. Another interpretation is that, instead of 
repeated geometric clearing of the line of sight to the central star and subsequent obscuration of 
this line by dust, the source is surrounded by dust that is being in situ reformed in the nearby 
circumstellar environment and then either destroyed or accelerated away though wind/outflow. In 
this scenario, an analogy to the deep occulting phenomenon of RCrB stars ( Claytori||2012| ) could be 



appropriate. The very large extinction variations, AAy up to 30, and the detected periodicity, may 
be a challenge for this hypothesis, however. The quasi-periodicity of the dust obscuration events in 
PTF lOnvg strongly suggest that the appropriate geometry is one of dense material along our line 
of sight to the optical/near-infrared continuum, that is in orbit about the central star, and located 
approximately 0.5-0.7 AU from it. 

The periodicity and the extinction-dominated nature of the 2011 and 2012 light curves do not, 
however extrapolate back to the 2010 (optical) data in a satisfactory manner, although the limited 
2009 data do seem to phase well with the derived period. During 2010 the source was much brighter 
than we have seen it over the subsequent two years, perhaps indicating that the absolute brightness 
peak was indeed driven by an accretion- related "outburst" event. The outburst hypothesis may be 
further supported by the lack of detection in previous infrared imaging surveys, notably 2MASS, 
but detection at not fainter than K s ~11 mag over the past three years. That there was brightening 
in the mid-infrared as well during 2010 suggests a bolometric luminosity increase, though the scale 
of the mid-infrared brightening could plausibly be due to the inferred extinction changes as well. 
Dramatic photometric and spectroscopic changes in young stars are often classified as outbursts 
and associated with large scale (one or more orders of magnitude) increases in the accretion rate of 
material from the disk on to the central star. Below we pursue the interpretation that the body of 



-57- 



evidence in hand for PTF lOnvg is indeed indicative of an outbursting young star, though mediated 
by semi-periodic dust obscuration events. 

The strong forbidden line emission is explained if the enhanced accretion episode of 2010 also 
drove an enhanced outflow, with the outflowing material reaching a shock-enabling surface approxi- 
mately one year post-outburst, thus causing the increase in the forbidden lines. However, given that 
we see the forbidden emission only when the source is in a faint state, and given the evidence for 
large variations in extinction towards the central continuum source, another interpretation is that it 
is only when the continuum is suppressed that we are able to see a much fainter, spatially distinct, 
but long-present jet region. This scenario requires an inclined disk to do the obscuring, but not too 
inclined given the >200 km/s line widths of the O I, S II, and Fe II lines. The situation would be 



similar to that advocated by White & Hillenbrand (2004): there is an orientation effect such that 



the spatially extended forbidden line emission is more readily observed in systems without a direct 
view of the central continuum source. This explanation was prefered to invoking differences, or in 
our case variations in, the mass outflow rate itself for explaining stronger forbidden line emission. 



We note that Appenzeller et al. ( 2005 ) also found enhanced line-to-continuum ratios for forbidden 



lines in known edge-on disk systems relative to those typically seen in accreting T Tauri stars. 
Another peculiarity of PTF lOnvg is the rare observation of optical TiO/VO molecular bands 



m emission. 



Hillenbrand et al. (2012) reported several additional sources besides PTF lOnvg that 
have been observed to exhibit similar, broad TiO/VO emission at optical wavelengths. In common 
between the sources is the Class I-type spectral energy distribution and spectroscopic evidence for 
rapid accretion and strong outflow. It is unclear at present whether an "outburst" phase is required 
for detection of TiO/VO emission. Future study of such TiO/VO sources may help elucidate the 
relative molecular abundances in circumstellar disk/envelope environments, and perhaps the mech- 
anisms for molecule formation and destruction. The TiO and VO species appear from the statistics 
to be not readily observable during more quiescent stages of disk evolution but notably they are 



present in planetary atmospheres and so must be present in the disks. Banzatti et al. (2012), for 
example, have recently discussed in detail for EX Lup the H2O and OH (increased emission during 
outburst), as well as several organic (decreased emission during outburst) molecules during quies- 
cence and in outburst. These authors discuss the photo-dissociative and photo-desorption effects 
of enhanced high energy radiation during an episode of enhanced accretion, and the ensuing broad 
implications for disk chemistry and structure evolution. 

In summary, a coherent physical model for PTF lOnvg would include the following. The optical 
TiO/VO and infrared CO emission comes from the inner disk region. The broad permitted atomic 
lines are formed in the accretion/outflow zone, with the narrower atomic lines part of the rest- 
frame accretion shock. The 2010 brightening was indeed an "outburst" event, driven by accretion, 
but perhaps accompanied by dust clearing along the line-of-sight. The photometric variability 
since 2010 is primarily due to large variations in extinction, which periodically obscure the inner 
disk. This causes the dense gas tracers (molecular TiO/VO, CO and permitted atomic lines) to 
dissappear. However, the low density tracers (molecular H2 and atomic forbidden lines) are coming 
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from an extended jet/outflow, which may include self-shocks. We can see this faint emission only 
when the central continuum source is obscured. The characteristic timescale of the extinction 
variations is few week occulations, with ~221 day periodicity. If these events are linked with an 
orbital timescales, it suggests that the occulting material is located at 0.5-0.7 AU. 



7.2. PTF lOnvg in the Context of other Large- Amplitude Variable Young Stars 

In this section we discuss the characteristics of PTF lOnvg relative to those of other young 
stellar objects that have exhibited outbursting behavior, finding some common elements with known 
classes, but also some unique aspects of PTF lOnvg. We advocate the position that our knowledge of 
the diversity of time- variable behavior in young stars is vastly incomplete, as is our understanding 
of the underlying physical mechanisms. Although FU Ori stars are certainly a distinguishable 
class once the candidates are suitably characterized, other variable categories such as EX Lup-type 
objects and V1647 Ori-type objects appear to contain a range of behaviors among their members. 
While all of these variable types are likely explained by time- variable accretion and related outflow, 
young star "eruptions" may manifest themselves in very different ways with the currently identified 
categories including somewhat diverse objects at present. 

The lightcurve of PTF lOnvg is dissimilar to those of other categories of young star outbursts. 
Beyond the initial rise, the light curve is very different from those of FU Ori stars, which remain 
at or near the elevated brightness for decades. The FU Ori stars are inferred to be undergoing 



significantly enhanced accretion at rates of roughly 10 -4 or 10 -5 M yr _1 (Zhu et al. 2010). The 
PTF lOnvg light curve is also dissimilar to the lightcurves of EX Lup-type objects, which have 
smaller amplitude rises and a clearly distinct quiescent vs outbursting state. The EX Lup star 
accretion rates are inferred to be lower than those advocated for FU stars; the unusually large 2008 
outburst of EX Lup itself was interpreted as an accretion enhancement from roughly W~ 9 Mqyt^ 1 
to roughly 10~ 7 Moyr" 1 (Sipos et al. 2009; Aspin et al. 2010; Juhasz et al. 2011). The lightcurve 
for this outburst event (see e.g. Figure 1 of Juhasz et al. 2012) indeed bears some similarity to the 
photometric behavior of PTF lOnvg, both in terms of the amplitude of the initial rapid rise and the 
multiple photometric peaks occuring over the next 1/2 to 3/4 of a year before the source faded to 
quiescent levels. As described here, PTF lOnvg also exhibits quasi-periodic brightness peaks at ~7 
month- long intervals. However, the outbursting phase with superposed large exinction variations 
has lasted almost 3 years now, longer than the events typically associated with EX Lup outbursts. 
It remains to be seen whether we are now witnessing the quiescent stage behavior, which would 
be the current stage of large amplitude extinction-driven periodic variability, or if a true quiescent 
stage will reveal itself in the future. 

The physical mechanism for the outbursting behavior in the several classes of large amplitude 
variable young stars is thought to be related to an inner disk instability. For FU Ori stars, thermal, 
gravitational, and magneto-rotational models have all been proposed (Bonnell &: Bastien 1992; 
Bell & Lin 1994; Clarke & Syer 1996; Kley & Lin 1999; Armitage et al. 2001; Vorobyov & Basu 



2005; Boley et al. 2006; Zhu et al. 2009b). For the EX Lup stars, the origin is postulated to be 
essentially at the inner disk edge itself, and the "instability" driven by a cyclic interaction between 
the differentially rotating accretion disk and the magnetosphere (roughly dipole in this region) which 
is rotating with the star. (Goodson h Winglee 1999 D'Angelo & Spruit 2010, 2012). Episodic 
accretion and outflow are both naturally explained as a consequence of the oscillations in magnetic 
field topology. The field structure alternates between an expanded magnetospheric structure which 
diffuses into the disk and leads eventually to field line opening and outflow. Following reconnection, 
a more compact magnetospheric structure results and allows encroachment by the accumulated 
inner disk material. Enhanced accretion on to the star results when co-rotation is crossed, after 
which the magnetosphere again expands. Alternately, the enhanced accretion could involve similar 



mechanisms but operating in a more stochastic rather than cyclic manner (e.g. Romanova et al. 



2005, 2008). In both scenarios, however, the predicted time scales are highly uncertain, ranging 



between several times the dynamical time scale and several orders of magnitude times the viscous 
time scale. Thus it is not yet known if the postulated inner disk oscillations can be associated 
with the interval between the young star outbursts that have been observed. At least qualitatively, 
there is some resemblance between the accretion rate evolution during the outbursts described in 
e.g. Figure 4 of D'Angelo & Spruit (2010) and the light curve behavior of PTF lOnvg. The same 
interpretative conclusion was reached for the 2008 outburst of EX Lup by Juhasz et al. (2012). 



Covey et al. (2011 ) made the analogy between the outburst behavior of PTF lOnvg and V1647 



Ori. This claim was before discovery of the large amplitude near-infared variations as consistent 



with extinction variations, and before discovery of the periodicity in the lightcurve. However, Aspin 



(2011) do claim that the near-infrared color variations of V1647 Ori are consistent with variable 
extinction (though they argue that those of PTF lOnvg are not consistent with this interpretation, 
contrary to the evidence presented here). Furthering the analogy, we have found a single outflow 
cavity in the long wavelength molecular line observations of PTF lOnvg which could also be argued 
for V1647 Ori based on the strong asymmetry in the scattered light nebula. 



Lorenzetti et al. (2012) have suggested based on the near-infrared color differences between 



the quiescent and the outbursting states of several EX Lup-type variables (including PTF lOnvg 
and V1647 Ori) that the outbursts are characterized by an extra thermal component that is not 
present during the quiescent stages. In the case of PTF lOnvg, the extra radiation during the 
outburst is estimated to have ~1400 K and to be located at 0.3 AU (based on a fitted blackbody 
model), after accounting for extinction variations of 4.6 mag between a low state and a high 
state. While our observations indicate that the observed photometric variations (Figure [6]) and 
also spectrophotometric variations (see later discussion of the de-reddened version of Figure 11 ) are 
roughly consistent with variations in extinction, there is some stochasticity to the color evolution 
which could indeed necessitate underlying changes in the source SED below the extinction variation 
which dominates. 



If the Lorenzetti et al. ( 2012 ) model is robust, the physical situation would be that of changes in 
thermally emitting material located in the inner region the disk. Although Lorenzetti et al. (2012) 
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interpret this as reprocessing of the radiation emitted on the stellar surface as the accretion shock 
becomes hotter due to the enhanced accretion rate, one can also draw analogy to the above theories 
invoking inner disk instabilities wherein the inner disk heating is due to viscous and dynamical 
process rather than radiative reprocessing. Notably, the ~0.3 AU scale derived in the Lorenzetti 



et al. (2012) blackbody model is not so different from the ~0.5-0.7 AU scale that was found above 
in our discussion of the semi-periodicity of the repeated photometric peaks. However, perhaps 
inconsistent between the postulation of an extra, cool emission component and the instability 
theory, is that there does not appear to be any lag between the infrared brightening and the optical 
brightening in PTF lOnvg. This is true for the initially recorded outburst in 2010 as well as all 
subsequent brightenings that we have observed from 2010 to 2012. Lorenzetti et al. (2012) do note 
that the blackbody model becomes less good of a fit as their EX Lup-type sources evolve beyond 
the initial photometric outburst towards quiescence. 



7.3. The Rate of PTF lOnvg-like events 



As detailed both here and in previous work (Covey et al. 2011 Aspin 2011), PTF lOnvg 



does not belong to either the FU Ori-like or the EX Lup-like classes of YSO outbursts. It may 
be an intermediate timescale and amplitude outburst class, perhaps analogous to V1647 Ori as 



discussed in Covey et al. (2011). Or it may be more a more extreme version of the extinction- 
driven phenomena seen in UX Ori or GM Cep. As discussed here, both physical phenomena 
probably play some role in the interpretation of PTF lOnvg. Estimates of the absolute rates of any 
such events remain highly uncertain. Nevertheless, we attempt a rough estimate of the relative rate 
of PTF lOnvg-type events from our ongoing photometric monitoring survey of the North America 
Nebula. History shows that astronomy is rich with examples of new phenomena being discovered 
as surveys push to explore undersampled regions of parameter space. 

To our knowledge, the P48 survey of the North America Nebula region is unique in both its 
depth and cadence. During the past ~3 yr the P48 has observed the same ~7.3 deg 2 field with 
roughly a nightly cadence while it was visible from Palomar mountain to a depth of ~20 mag. 



In this time there has been one FU Ori-like outburst (HBC 722/PTF lOqpf; Semkov et al. 2010 



Miller et al. 2011) and one PTF lOnvg-like event, suggesting that the rate of the two phenomena 
is similar to within the accuracy afforded by low number statistics. Overall, with ~12 FU Ori-like 



eruptions observed to date (Reipurth k, Aspin 2010) and only one PTF lOnvg-like event it is also 
possible that PTF lOnvg-like events are far more rare and only ~one-tenth as common as FU Ori 
outbursts. With ~10 FU Ori-like outbursts expected for each low mass star (Hartmann & Kenyon 



1996), even this lower limit would suggest ~1 PTF lOnvg-like outburst per low mass star. Likewise, 



in the interpretation that this event corresponds to the initial cloud clearing and source revelation, 
we would expect ~1 event per star. 



As future surveys continue to explore star forming regions with high cadence observations over 
long time baselines, the importance of PTF lOnvg-like events in our global understanding of the 
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low mass star formation process will be elucidated. They may become common in surveys adopting 
a similar strategy, or they may be revealed as rare events. 

8. Summary and Conclusions 

1. Continued photometric monitoring of PTF lOnvg by the PTF 48" and PAIRITEL telescopes 
throughout 2011 and 2012 detected magnitude variations of AR ~ 6 mag and AK ~ 3 mag, and 
color variations of AJ-K ~ 3 mag. The time series behavior in color-color and color-magnitude 
diagrams indicates that these photometric variations are consistent with those expected due to 
changes in line-of-sight extinction amounting to AAy > 30 mag. 

2. Over the long term, attributing PTF lOnvg's non-detection in 2MASS to similar extinction 
changes suggests the source has historically experienced a range of AAy ~ 70 mag. Alternately, 
the non-detection in 2MASS may reflect a true outburst scenario for the 2010 brightening and 
subsequent photometric evolution. 

3. Time scales of several weeks can be associated with dramatic flare- like photometric changes 
and ~7 months with the repeated maxima observed in the time series over the past several years. 
While the initial (2010) maxima were rounded/smooth, later maxima were characterized by nar- 
row/sharp photometric peaks. 

4. The Lomb-Scargle periodogram calculated from JHK S light curves features a prominent, 
statistically significant peak at ~221 days. Continued monitoring will be essential to confirm this 
signal, however, as the existing light curves span <1000 days in total, and include several 100-1- 
day gaps. If this signal persists, it suggests under the assumption that the period is the Keplerian 
orbital timescale, that the dust obscuration governing PTF lOnvg's photometry is located ~0.5-0.7 
AU from the central star. 

5. WISE detected PTF lOnvg over multiple epochs in 2010. Comparison of these mid-IR 
detections with prior detections at these wavelengths by IRAS, Spitzer, and MSX indicate that 
PTF lOnvg brightened in the mid-infrared by factors of a few during its summer 2010 brightening. 

6. CARMA observations reveal an extended millimeter continuum source as well as spatially 
unresolved emission in both 12CO and 13CO centered on PTF lOnvg's near-infrared position and 
systemic velocity. Spatially extended, redshifted 12CO emission is also detected to the south of 
PTF lOnvg, which we interpret as revealing the presence of an outflow cavity. Over the wider field, 
the molecular emission is well-aligned with the Ha emission arc. 

7. The bright-state optical and near-infrared spectrum of PTF lOnvg is that of a "continuum- 
plus-emission" object (with significant molecular emission contributions to the continuum) whereas 
the faint-state spectrum is that of a Herbig-Haro object. Both neutral and singly ionized line 
species probing a range of densities and temperatures are observed. The lines that characterize 
PTF lOnvg's faint-state spectrum are present at all epochs, however, suggesting that the spectral 
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evolution is primarily due to the supression of the bright state spectral features during epochs of 
enhanced extinction. 

8. Spectral monitoring reveals that the strong wind signatures detected in the 2010 outburst 
persisted through 2011 and 2012. The depth and terminal velocities of optical wind absorption 
features (i.e., the Na I D and K I 7665/7699 doublets, the O I 7773 triplet) appeared consistent 
across all epochs where the optical continuum was sufficiently strong to enable these features to be 
measured. Moderate resolution near-infrared spectra of PTF lOnvg indicate that the He I 10830 
line increased in both depth and equivalent width by nearly a factor of two from 2010 through 
2011. High- resolution He I 10830 line profiles are consistent with this long term trend. A strong 
redshifted emission component to this line likely also traces the outflow rate of the inner wind. 

9. The evolution of the strengths and profiles of permitted atomic emission lines reveal the 
presence of time variability in PTF lOnvg's line-of-sight accretion and outflow activity. The source 
exhibits optically thick Ca II triplet emission with line ratios that are unusual even for young stars. 
The Ha emission also demonstrates significant variations in both overall strength and kinematic 
profile, with the line weakening and shifting blueward as the source fades. Strong redshifted Ha 
emission disappeared in June 2011 as PTF lOnvg approached its photometric minimum, while the 
blueshifted Ha emission remained distinctly visible. These observations suggest that PTF lOnvg's 
Ca II and H I lines include contributions from both accretion and outflow activity. While the 
accretion component may dominate during the brightest states, during fainter states only the 
outflow signature remains. Notably absent from the spectra at any epoch is the He I A 5876 line, 
a high temperature (T > 15,000 K) tracer seen near-ubiquitously in classical T Tauri stars. 

10. Forbidden line emission measuring ~5000-20,0000 K gas reveals the presence of a likely 
Herbig Haro jet with a wide range of densities, n ~ 10 3 — 10 7 cm -3 . Some line ratios such as [O I] 
6300 A : 6363 A and [Fe II] 7155 A : 7172 A show extreme and unusually time-dependent ratios 
that are not easily explained in the standard optically thin assumption for forbidden line emission. 
Instead, these lines may be partially optically thick. 

10. Spatially extended blueshifted emission is visible in several forbidden line species as well as 
Ha in two-dimensional HIRES images. These position-velocity diagrams provide further evidence 
that PTF lOnvg is driving a jet with a line-of-sight velocity of about —200 km/s and exciting a 
strong shock region ~4" projected distance from the source (2000AU at 520 pc). 

11. The moderate resolution near-infrared spectra obtained in June/ July 2011, when PTF lOnvg 
was near its faintest state, provide an opportunity to study the properties of the jet and extended 
outflow. PTF lOnvg's H2 line strengths strongly resemble theoretical predictions for C-type shocks 
and pure X-ray photo-excitation, as well as empirical observations of H2 emission from RW Aur 
and Knot B of the HH54 protostellar jet; the same is true for the near-infrared [Fe II] emission. 

12. From near-infrared adaptive optics observations we can exclude the presence of stellar 
companions within several hundred AU of PTF lOnvg. 
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Further study of this enigmatic source is warranted. 

Facilities: Palomar 48", PAIRITEL, Lick 3m, Palomar 200" (two instruments), APO, Keck 
(four instruments), IRTF, CARMA, DSS, KPNO 0.9m, USNO, 2MASS, UKIDSS, Spitzer 
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Table I. P48 iZ-Band Photometry of PTFlOnvg 1 



Epoch (MJD) K, PTF (mag) err 



55056.221 17.831 0.019 

55056.295 17.772 0.020 

55059.262 17.328 0.017 

55059.357 17.380 0.017 

55061.316 17.520 0.023 



lr The table will appear in its entirety in the 
electronic edition of the journal. 



Table 2. PAIRITEL JHK-Band Photometry of PTFlOnvg 1 



Epoch (MJD) 


J mag 


,/ err 


H mag 


H err 


K mag 


K err 


55387.376 


12.229 


0.025 


10.232 


0.030 


8.456 s 


0.060 s 


55392.294 


12.231 


0.026 


10.231 


0.030 


8.475 s 


0.063 s 


55450.191 


10.896 


0.025 


9.414 s 


0.031 s 


8.029 s 


0.060 s 


55463.151 


11.538 


0.025 


9.801 


0.030 


8.303 s 


0.060 s 


55466.126 


11.585 


0.025 


9.966 


0.030 


8.347 s 


0.060 s 



"'"The table will appear in its entirety in the electronic edition of the journal. 

" s Photomctry measured on short read PAIRITEL images (see text). 

"No reliable photometric measurements arc available due to large thermal back- 
grounds (sec text). 
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Table 5. PTF lOnvg [Fe II] Line Fluxes 



Lino ID 


PTF lOnvg 


G06 HH54 


1.257 Mm 


4.5 


8.7 


1.644 |jm 


5.06 


8.2 


1.257 / 1.644 


0.89 


1.06 



Table 6. PTF lOnvg H 2 Line Fluxes 



Line ID 


PTF lOnvg 


B08 RW Aur (spatially extended) 


S95 C Shock 


S95 J Shock 


G06 HH54 


N07 X-ray (10 Mm) 


N07 UV (10 Mm) 


1.958 fim [1-0 S(3)] 


1.8 




2.68 


0.0341 


15.4 


1.326 


8.096 


2.034 fim [1-0 S(2)j 


0.75 


0.25 


1.06 


0.011 


4.2 


0.533 


3.857 


2.073 fim [2-1 S(3)j 


0.192 


0.12 


0.139 


0.00868 


1.63 






2.123 pm [1-0 S(l)l 


2.045 


1.0 


3.04 


0.0266 


9.59 


1.565 


12.968 


2.154 Mm [2-1 S(2)j 


0.041 


<0.09 


0.0539 


0.0027 


0.51 






2.222 Mm [1-0 S(0)j 


0.476 


0.17 


0.679 


0.00526 


2.58 


0.357 


3.256 


2.248 Mm [2-1 S(l)j 


0.206 


0.11 


0.151 


0.00635 


1.00 


0.09 


0.242 


2.407 Mm [1-0 Q(l)] 


1.63 


0.76 


2.34 


0.0167 


10.9 


1.27 


12.299 


2.413 Mm [1-0 Q(2)j 


0.557 


0.25 


0.749 


0.0058 


4.5 


0.395 


3.604 


2.248 / 2.123 


0.10 


0.11 


0.05 


0.24 


0.10 


0.06 


0.02 


2.034 / 2.222 


1.58 


1.52 


1.56 


2.09 


1.63 


1.49 


1.18 


2.248 / 2.073 


1.07 


1.09 


1.09 


0.73 


0.61 






2.123 / 2.407 


1.25 


1.38 


1.30 


1.59 


0.88 


1.23 


1.05 
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